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ABSTRACT 


The  Applied  Physies  Laboratory  of  the  University  of  Washington  (APL-UW)  aequired 
on  loan  an  experimental  deviee  known  as  the  multiport  transdueer.  APL-UW  developed,  in 
turn,  a  eomplete  transmitter  system  that  integrates  this  transdueer,  eapable  of  wideband  op¬ 
eration  (roughly  180-350  Hz)  from  near-surfaee  depths  to  depths  greater  than  1000  m.  The 
system’s  eleetrieal  eomponents  inelude  an  autotransformer  tuner,  a  battery  power  module, 
and  a  fibre  optie  telemetry  interfaee;  meehanieal  eomponents  inelude  a  steel  supporting 
strueture  and  a  pressure-eompensated  tuner  housing;  an  additional  aeoustieal  eomponent 
is  a  monitor  hydrophone  in  a  vibration  isolation  mount;  and  a  signal  eomponent  involves 
a  lumped  parameter  SPICE  eireuit  model  approximation  of  the  entire  end-to-end  system, 
an  assoeiated  C-f-f  applieation  to  prediet  the  time-domain  aeoustie  far  field  from  a  standard 
time-domain  waveform  input  file,  and  a  pre-equalization  filter.  The  multiport  system  was  a 
key  element  in  a  2009  at-sea  oeean  aeousties  experiment  loeated  in  the  Philippine  Sea  and 
provided  many  hours  of  high-quality  pulsed  transmissions  to  a  nearby  vertieal  line  array  of 
hydrophones. 
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EXECUTIVE  SUMMARY 


The  Applied  Physies  Laboratory  of  the  University  of  Washington  (APL-UW)  aequired 
on  loan  an  experimental  deviee  known  as  the  multiport  transdueer.  This  deviee  is  a  double 
ported  free-flooded  tube  resonator,  with  resonanees  at  about  210  Hz  and  320  Hz.  APL-UW 
developed,  in  turn,  a  eomplete  transmitter  system  that  integrates  this  transdueer,  whieh  is 
deseribed  in  this  report. 

The  inventor  of  the  multiport  transdueer  (Dr.  Jaek  Butler,  ImageAeousties,  Ine.)  main¬ 
tains  a  finite-element  model  of  the  transdueer,  but  for  the  purposes  of  the  effort  deseribed 
here,  a  simpler  lumped  parameter  equivalent  eireuit  was  developed  that  models  the  transfer 
of  input  voltage  to  “output”  fluid  veloeity  at  the  openings  of  the  tubes.  The  system  parame¬ 
ters  were  determined  based  on  eleetrieal  and  aeoustieal  eharaeteristies  measured  during  an 
engineering  test  eondueted  by  APL-UW  in  Lake  Washington  in  2006. 

Transdueers  are  usually  mated  to  matehing  tuners,  and  several  tuners  were  eonsidered. 
In  addition  to  tuning  out  the  impedanee  mismateh  at  various  frequeneies,  the  tuners  also 
served  to  step  up  the  voltage  into  the  transdueer.  This  reduees  the  neeessary  voltage  rating 
on  the  drive  eable. 

All  the  signal  eonditioning  eomponents  from  the  digital-to-analog  eonverter  to  the 
transdueer  were  implemented  as  subeireuits  in  a  SPICE  model  (Simulation  Program  with 
Integrated  Cireuit  Emphasis;  [I]),  and  ineorporated  into  a  C-f-f  applieation  eode  that  takes 
an  input  file  eontaining  a  time-domain  waveform  (eneoded  into  12-bit  digital-to-analog 
words)  and  eomputes  the  time-domain  far  field  pressure  waveform.  Various  uses  of  this 
eode  are  demonstrated:  the  RMS  sealing  differenee  between  broadside  souree  level  and  an 
endfire  monitor  sensor  level,  the  wide  band  beam  pattern,  and  aetual  souree  level  ealeula- 
tions. 

The  transfer  funetion  into  far  field  broadside  pressure  is  signifieantly  not  flat  but  double- 
peaked,  so  there  is  a  need  to  equalize  the  drive  signal  to  mitigate  the  phase  and  amplitude 
distortion  produeed  by  the  transdueer.  A  pre-equalization  filter  was  developed,  and  sim¬ 
ulations  demonstrate  the  effleaey  of  pre-equalization  (based  on  the  SPICE  model  transfer 
funetion)  on  pulse-eompressed  m-sequenees. 

The  multiport  system  was  a  key  element  in  a  2009  at-sea  oeean  aeousties  experiment 
in  the  Philippine  Sea  and  provided  many  hours  of  high-quality  pulsed  transmissions  to  a 
nearby  vertieal  line  array  of  hydrophones.  A  listing  of  these  transmissions  is  provided  in 
Appendix  C. 
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1  Introduction 


In  2005  the  Applied  Physics  Laboratory  of  the  University  of  Washington  (APL-UW)  ac¬ 
quired  on  loan  an  experimental  device  known  as  the  “multiport”  transducer  from  the  Naval 
Undersea  Warfare  Center  Division  Newport.  This  transducer  was  originally  designed  and 
constructed  by  ImageAcoustics,  Inc.  [Cohasset,  MA]  and  Massa  Products,  Inc.  [Hingham, 
MA]  under  a  Small  Business  Innovative  Research  (SBIR)  grant  to  ImageAcoustics,  Inc. 
from  the  Naval  Undersea  Warfare  Center  Division  Newport. 

APL-UW  was  subsequently  funded  by  Code  3210A  of  the  Office  of  Naval  Research 
(ONR)  to  develop  a  deep-water  ship-suspended  system  based  around  the  MP200/TR1446 
for  operations  in  the  ONR-sponsored  Philippine  Sea  experiments  of  2009  and  2010.  This 
report  summarizes  that  development  effort.  Section  2  describes  the  development  effort  for 
a  simple  circuit  model  suitable  for  routine  system  performance  calculations.  Some  of  these 
results  have  been  disseminated  in  previous  notes;  the  calculations  regarding  the  transfor¬ 
mation  ratio  are  new.  Chapter  3  describes  the  efforts  made  to  design  and  build  a  tuner  for 
the  transducer.  Section  4  describes  the  entire  multiport  system  as  built  and  deployed  by 
APL-UW.  Section  5  describes  an  entire  system  model  of  the  multiport  system,  including 
the  top-side  electronics  and  the  electro-opto-mechanical  suspension  cable,  and  the  prop¬ 
agator  to  predict  the  far  field  time-domain  radiated  acoustic  pressure.  This  model  was 
incorporated  for  ease  of  use  into  a  simulator  code  that  takes  as  input  a  time-domain  “drive” 
waveform  in  convenient  standard  .  wav  file  format.  Section  6  discusses  the  implementation 
of  two  different  equalizers,  necessary  to  compensate  for  the  considerably  non-flat  response 
of  the  transducer  itself.  Finally,  section  7  provides  a  review  of  the  system  in  use  during  an 
at-sea  deployment  on  an  acoustics  experiment  in  the  Philippine  Sea  in  2009. 
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2  The  Multiport  Transducer 

2.1  Introduction 


The  MP200/TR1446  is  an  experimental  aeoustie  transdueer  invented  and  designed  by  Image- 
Aeousties,  Ine.  [Cohasset,  MA]  and  built  by  Massa  Produets,  Ine.  [Hingham,  MA;  Fig.  1]. 
The  theory  of  operation  is  fully  deseribed  in  the  original  patent  [2] .  In  brief,  it  is  a  dou¬ 
ble  ported  free-flooded  tube  resonator,  with  resonanees  at  about  210  Hz  and  320  Hz.  This 
deviee  does  not  require  pressurization,  and  ean  be  driven  to  radiate  at  least  195  dB  (nar¬ 
rowband)  at  any  depth. 

ImageAeousties,  Ine.,  maintains  a  eomplete  model  for  this  transdueer,  ineluding  a 
finite-element  eomponent  to  model  the  aeoustie  field  inside  the  tubes.  For  the  purposes 
of  the  development  of  the  APL-UW  system,  it  was  not  neeessary  to  utilize  their  complete 
model.  A  simpler  approximate  model  was  developed  at  APL-UW  and  is  presented  in  this 
section. 


Figure  1:  The  MP200/TR1446.  The  orientation  shown  has  the  axis  of  the  two  tubes  hori¬ 
zontal.  From  Refs.  [3],  [4]. 
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i  'rad 


Figure  2:  Simple  equivalent  electrical  circuit  for  the  Multiport  Source.  The  “currents”  i,  io 
and  iRAD  are  defined  in  section  2.5. 

2.2  Equivalent  Electrical  Circuit  —  Theory 

Characterization  of  the  equivalent  electrical  circuit  parameters  follows  the  procedures  of 
Wilson  [5],  Stansfield  [6],  and  Sherman  and  Butler  [7].  Given  a  power  amplifier  drive 
voltage  of  v{t)  with  Fourier  transform  V (/)  and  a  drive  current  i(t)  with  Fourier  transform 
/(/),  the  admittance  is 

Y(f)^^^G(f)+jB(f)  (1) 

where  G{f  )  is  the  conductance  and  B{f  )  is  the  susceptance,  both  in  units  of  siemens,  and 
j  =  y/^.  The  curves  G{f),  B{f)  and  the  admittance  locus  Y  (/)  can  be  used  to  determine 
various  “raw”  parameters,  as  described  below,  which  are  then  converted  to  the  elements  of 
a  simple  equivalent  circuit. 

A  simple  but  adequate  equivalent  circuit  diagram  (Fig.  2)  represents  the  blocked  capac¬ 
itance  of  the  device  with  a  shunt  capacitor,  and  the  two  tubes  of  the  device  with  two  parallel 
series  RLC  branches.  There  are  two  distinct  resonances  in  the  device,  each  represented  by 
an  RLC  branch  (Fig.  2),  and  these  are  far  enough  apart  in  frequency  space  that  a  sepa¬ 
rate  analysis  is  pursued  for  each.  The  raw  parameters  and  their  computation  are  described 
below: 


Low-Frequency  Capacitance.  The  effective  input  admittance  at  low  frequency  goes  as 
B{uj)/uj  as  u  goes  to  zero  and  is  represented  by  the  capacitance  Grr. 

Motional  Resonance.  The  motional  resonance  frequency  fs  is  determined  by  locating  the 
peak  in  the  conductance  G{f  ). 
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Mechanical  Q.  The  mechanical  Q,  denoted  Qm,  is  determined  by  the  equation 

Qm  =  -T^  (2) 

h  ~  ji 

where  /o  is  a  resonance  frequency,  and  {fi,  f 2)  are  the  frequencies  where  the  conduc¬ 
tance  drops  to  half  its  maximum  value  at  the  resonance  peak.  These  half-amplitude 
frequencies  were  calculated  using  linear  (inverse)  interpolation  about  G  =  0.5. 

Parallel  Resonance.  The  parallel  resonance  frequency  fp  is  determined  by  finding  the 
frequency  fp  >  fs  where  the  Y (/)  locus  meets  the  line  from  the  origin  to  the  point 
{G{fs),  B{fs))  in  the  admittance  plane. 


Once  these  parameters  have  been  determined,  the  following  circuit  parameters  can  be 
calculated: 


Coupling  Coefficient.  The  (squared)  coupling  coefficient  is 

=  I- $  (3) 

fp 

Mechanical  Capacitance.  The  effective  resonator  mechanical  capacitance,  denoted  by  G, 
is  roughly 

C  =  ^Glp  (4) 

Resonance  Inductance.  The  effective  resonator  mechanical  inductance  L  is  then  related 
to  G  via  the  resonance  frequency: 


L  = 


ulG 


(5) 


Blocked  Capacitance.  The  blocked  admittance  is  assumed  here  to  be  purely  capacitative, 
and  the  capacitance  is 

Go  =  Glf  —  G  (6) 


Resonance  Resistance.  The  resistance  at  resonance  R  is  given  roughly  by 

R=  ^ 


^^sCQm 


(7) 
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2.3  2006  Lake  Washington  Test  Configuration 

APL-UW  tested  the  MP200/TR1446  over  the  side  of  the  R/V  Henderson  for  several  days  in 
April  2006  in  Lake  Washington.  The  transducer  was  driven  by  an  80486  computer  running 
MS-DOS  and  the  program  atocsam.  The  drive  signal  from  the  computer  was  routed  to 
a  Ling  power  amplifier  and  via  the  suspension  cable  to  the  transducer.  A  custom  circuit 
internal  to  the  Ling  monitors  the  “top-of-the-cable”  drive  voltage  and  current. 

An  ITC  6050C  hydrophone  was  deployed  over  the  side  of  the  R/V  Henderson  to  moni¬ 
tor  any  acoustic  emissions  from  the  transducer.  The  hydrophone  signal  was  routed  into  the 
“Tennelec”  amplifier  box  (approximately  no  gain)  through  a  Krohn-Hite  filter  (dialed  to 
pass  100  Hz  to  400  Hz  at  0  dB),  through  an  HP  4436A  attenuator  (set  to  0  dB  attenuation), 
and  then  into  the  “bare”  signal  input  on  the  card  cage  chassis. 

The  input  signal  was  generated  via  the  application  seqsam  and  consisted  of  an  m- 
sequence  with  octal  law  2033.  The  modulation  angle  was  88.2092°,  2  cycles  per  digit, 
1023  digits.  The  carrier  frequency  was  250  Hz.  As  is  customary  with  ATOC  (acoustic 
thermometry  of  ocean  climate)  signals,  the  waveform  was  oversampled  by  a  factor  of  32, 
so  the  sampling  rate  of  the  playback  digital-analog  board  was  8  kHz. 

The  program  atocsam  also  performs  simultaneous  analog-to-digital  while  stream¬ 
ing  out  through  the  digital-to-analog  (D/A)  converter  subsystem.  In  normal  operation, 
five  channels  are  captured:  the  1  pulse/second  line,  power  amplifier  drive  voltage  moni¬ 
tor,  power  amplifier  drive  current  monitor,  “bare”  acoustic  input,  and  IRIG-M.  Due  to  the 
high  sampling  rate,  it  was  not  possible  to  capture  all  five  channels,  and  perform  D/A,  with 
the  computer  power  available.  Therefore,  only  three  channels  were  captured:  the  voltage 
monitor,  the  current  monitor,  and  the  bare  acoustic  line. 

Each  test  measurement  was  120  s  in  duration.  The  programmable  anti-alias  filters  on 
the  filter/gain  card  were  set  to  384  Hz  automatically  by  atocsam. 


2.4  Equivalent  Circuit  Characterization  —  Results 

With  the  Ling  power  amplifier,  the  voltage  monitor  circuit  divides  the  power  amplifier  drive 
voltage  by  a  factor  of  1000,  and  the  current  monitor  circuit  converts  the  power  amplifier 
drive  current  to  a  voltage  by  200  mV/A. 

In  practice,  V (/)  and  1(f)  are  computed  using  discrete  Fourier  transforms.  In  regions 
outside  the  main  band  of  the  signal,  these  terms  are  largely  noise,  and  therefore  only  the 
admittance  over  the  main  frequency  band  /  G  [/min,  /max]  need  be  considered. 
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Figure  3:  B{uj)/uj  versus  frequency,  2006  Lake  Washington  test.  File  mptestOl .  dat. 

In  a  120-s  transmission,  there  are  14  whole  waveforms.  Therefore,  to  mitigate  additive 
system  noise  and  increase  the  reliability  of  the  admittance  estimates,  the  14  m-sequence 
waveforms  were  coherently  averaged  prior  to  the  discrete  Fourier  transform  step. 


Low  Frequency  Capacitance.  The  effective  input  admittance  at  low  frequency  goes  as 
B{uj)/uj  as  uj  goes  to  zero  and  is  represented  by  the  capacitance  Clf-  In  practice, 
the  admittance  loci  were  only  calculated  down  to  100  Hz  (because  signal  to  noise 
issues  progressively  contaminated  the  result  at  lower  frequencies).  Furthermore,  to 
moderate  the  effect  of  noise,  the  lowest  five  values  of  B{uj)/uj  were  averaged.  As  an 
example,  see  Fig.  3. 

Motional  Resonance.  The  motional  resonance  frequency  fs  is  determined  by  locating  the 
peak  in  the  conductance  G{f).  (The  function  was  represented  by  a  tab  listing  sam¬ 
pled  at  122  mHz:  here  the  frequency  of  the  nearest  frequency  grid  point  was  used.) 

Mechanical  Q.  The  mechanical  Q  is  determined  by  Eq.  2,  reproduced  for  convenience 
below: 

p  fo 

VM  —  7 - J- 

h  —  h 

where  fo  is  a  resonance  frequency,  and  {fi,  f 2)  are  the  frequencies  where  the  conduc¬ 
tance  drops  to  half  its  maximum  value  at  the  resonance  peak.  These  half-amplitude 
frequencies  were  calculated  using  linear  (inverse)  interpolation  about  G  =  0.5. 
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Figure  4:  Fourier  transforms  of  drive  signals,  2006  Lake  Washington  test.  Top:  drive  signal 
from  file  mpO  0 1 .  sam.  Middle:  power  amplifier  drive  voltage  from  file  mptest  1 5  .  sam. 
Bottom:  P.A.  drive  current  from  file  mptest  15  .  sam.  The  middle  and  bottom  panels 
utilize  the  first  m-sequence  in  the  file. 


The  Fourier  transforms  of  the  drive  signal  and  the  voltage  and  current  monitor  channels 
for  one  m-sequence  are  shown  in  Fig.  4.  The  MP200/TR1446  has  two  resonances  at  ap¬ 
proximately  200  and  300  Hz.  Note  that  there  is  a  smoothing  filter  after  the  D/A  prior  to  the 
Ling  to  prevent  high-frequency  energy  from  entering  the  amplifier,  and  this  can  be  seen  in 
Fig.  4  causing  an  attenuation  of  the  upper  half  of  the  main  lobe  and  the  upper  sidelobe  in 
the  drive  voltage  transform.  (This  should  have  no  effect  on  the  admittance  calculation  be¬ 
cause  it  causes  an  identical  attenuation  in  both  the  drive  voltage  and  current,  and  therefore 
cancels  in  their  ratio.) 

The  computation  of  these  parameters  from  an  input  admittance  file  was  automated  in 
a  Matlab  function,  and  applied  to  all  the  admittance  files  from  the  in-water  2006  Lake 
Washington  test  (with  the  exception  of  files  corresponding  to  unusable  transmissions.) 
Twenty-one  transmissions  were  conducted.  A  summary  of  each  collection  is  provided  in 
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Figure  5:  Admittance  locus,  SPICE  simulation. 


Appendix  A.  Summaries  of  the  measured  raw  and  computed  values  themselves  are  shown 
in  Tables  1  and  2. 

The  admittance  locus  produced  by  a  SPICE  simulation  is  shown  in  Eig.  5.  It  bears  good 
resemblance  to  the  measured  admittance  loci  (see  Eig.  65  in  Appendix  A,  for  example). 


2.5  Electro-Acoustic  Transformation  Ratio 

To  model  the  transducer  from  input  voltage  to  output  radiated  acoustic  field,  a  further  re¬ 
finement  to  the  equivalent  circuit  model  is  needed.  The  standard  approach  is  to  utilize  two 
transformation  elements,  one  to  transform  between  the  electrical  and  mechanical  stages, 
and  the  second  to  transform  between  the  mechanical  and  acoustic  stages.  The  approach 
taken  here  is  much  simpler:  a  single  transformer  is  used  to  convert  from  the  electrical  stage 
to  the  acoustic  field.  This  simplification  is  motivated  by  the  lack  (in  this  model)  of  internal 
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S  :  1 


Figure  6:  Electro-acoustic  transformation  element  for  the  impedance  analogy,  converting 
from  voltage  and  current  in  the  electrical  stage  to  force  and  velocity  in  the  acoustic  stage. 

— ^55855^ — ^ - • 

cfo  Ho  R’lo  ““ 


Figure  7 :  The  modified  equivalent  circuit  involving  the  transformation  element. 


mechanical  and  acoustic  variables,  and  the  need  to  only  match  input  voltage  levels  to  mea¬ 
sured  output  acoustic  levels.  This  “electro-acoustic  transformer”  is  shown  in  Fig.  6.  The 
modeling  goal  will  be  to  determine  the  transformation  ratio  S,  sometimes  called  the  turns 
ratio. 

First  consider  the  lumped  impedance  analogy  circuit  (Fig.  7),  which  now  includes  the 
electro-acoustic  transformer.  At  the  resonance  frequency  of  one  of  the  branches,  the  ve¬ 
locity  u  in  the  transformer  secondary  will  be  dominated  by  the  velocity  “flowing”  through 
that  branch.  The  problem  is  separated  into  two  distinct  problems,  one  at  each  resonance 
frequency.  Using  the  standard  rules  for  transformers,  the  current  and  velocity  are  related 
according  to 

u  =  Sirad  =  OiSi  (8) 

where  a  is  the  proportion  of  i  that  does  not  flow  through  the  blocked  capacitance  Cq.  This 
proportionality  can  be  determined,  e.g.,  from  SPICE. 

Next,  consider  the  radiated  acoustic  field  from  one  tube.  The  acoustic  wavelength  at 
resonance  is  much  greater  than  the  radius  of  the  opening  of  that  tube.  The  far  field  can  be 
approximated  to  first  order  as  if  it  originated  from  a  “simple  source”  monopole  located  at 
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that  end  of  that  tube.  This  is  true  for  both  tubes.  Thus,  let 


PM(r) 


47rr 


(9) 


be  the  radiated  pressure  due  to  a  monopole,  where  |r|  =  r  and  i  =  The  source 

strength  IIm  is  related  to  the  fluid  velocity  of  the  monopole  via  [8] 


YIm  =  pckU 


(10) 


where  k  =  u/c  and  U  is  the  volume  velocity  of  the  source.  For  a  simple  source  with 
radius  a,  the  volume  velocity  is  Ana^Ur,  where  Ur  is  the  radial  velocity.  A  simple  spherical 
source  is  not  conceptually  accurate  for  the  MP200/TR1446;  a  better  model  might  use  a 
piston,  with  volume  velocity  U  =  na^Uz,  where  Uz  is  now  the  velocity  of  the  fluid  at  the 
opening  of  the  tube,  assumed  constant  over  the  opening,  normal  to  the  opening.  The  fluid 
velocity  Ur  or  Uz  is  provided  by  the  SPICE  model:  the  choice  of  simple  source  or  piston 
volume  velocity  is  merely  a  matter  of  a  constant  of  proportionality,  but  the  choice  must 
be  consistent  with  the  propagator  used  to  calculate  the  far  field  (section  4).  To  maintain 
fidelity  with  the  physical  model. 


Hm  =  Tia^pckuz  (11) 

The  two  ends  of  the  tube  form  an  ‘acoustic  doublet,’  shown  schematically  in  Fig.  8. 
The  total  acoustic  field  at  the  field  point  is  the  superposition  of  the  contributions  from  both 
sources.  For  measurements  in  the  equatorial  plane  of  the  tube  (i.e.,  r  =  [xh,  0,  0]),  the 
two  contributions  are  in-phase.  Thus,  if  xh  ^  L  (where  L  is  the  length  of  the  tube),  the 
measured  acoustic  pressure  is 

\Ph{xhW  -  Mpm{xh)\‘^  (12) 


Hence,  using  Eqs.  9,  11,  and  12, 


4a;^ 


Uz  = 


H 


a?{pcY{kaY 


\Ph{xh)\^ 


(13) 


Equating  Eqs.  8  and  13,  the  far  field  pressure  can  be  written  as 


\Ph{xh)?  = 


(14) 


i.e.,  the  pressure  is  proportional  to  the  drive  current.  The  constant  of  proportionality,  de¬ 
noted  here  by  B,  consists  of  several  factors: 


B^ 


a?{pcy(ka)‘^a^S‘^ 

Axl 


(15) 
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Z 


Figure  8:  Acoustic  doublet.  A  resonating  tube  or  length  L  oriented  vertically  is  represented 
by  two  in-phase  equi-amplitude  monopoles  at  (0,0,+L/2)  and  (0,0,— L/2).  The  field 
measurement  is  made  at  location  {xh,  0,0),  where  xh  ^  L. 


If  S  were  known  from  theory,  B  could  be  computed,  and  so  too  the  far  field  pressure. 
This  would  be  a  prediction.  The  approach  here  is  to  utilize  empirical  measurements  of  the 
pressure  and  drive  current  to  estimate  S:  Eq.  14  gives 


^2  ^  \Ph{xh)\‘^ 


(16) 


In  principle,  this  is  adequate  to  compute  the  pressure  modulus  at  one  location  ((x^,  0,  0)) 
and  at  one  frequency.  This  ratio  is  likely  to  be  frequency-dependent,  and  tabulating  B{f)  as 
a  function  of  frequency  is  inconvenient.  Rather,  assuming  that  S  is  frequency-insensitive, 
we  therefore  seek  to  determine  it  explicitly.  Thus, 

o2  \Ph{xh)\^ 

a^a^{pcy{kaY 


To  proceed  further,  it  is  required  to  estimate  the  drive  current  divider  proportionality 
factor  a,  the  range  xh,  ka,  the  pressure  magnitude  at  the  field  point  \p{xh)\'^,  and  the 
corresponding  drive  current  magnitude  |. 


2.5.1  The  Parameter  ka 

The  dimensionless  parameter  ka,  where  a  is  the  radius  of  either  tube,  is  shown  in  Table  3  for 
both  tubes.  The  full  radius  of  the  outer  tube  is  used  for  the  corresponding  a  for  simplicity, 
although  the  radiating  region  is  annular,  and  not  tto^.  ka  is  less  than  unity  in  both  cases. 
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tube 

diameter 

resonance 

length 

[m] 

[Hz] 

ka 

long 

0.26 

212 

0.31 

short 

0.39 

320 

0.59 

Table  3:  Values  of  ka  for  both  tubes. 


file 

peak  bin 

delay  [ms] 

range  [m] 

mptest 04 . dat 

172 

21.5 

30.7 

mptest 05 . dat 

173 

21.6 

30.9 

mptest 0  6 . dat 

172 

21.5 

30.7 

Table  4:  Slant  ranges  from  transducer  to  monitor  hydrophone,  2006  Lake  Washington  test. 
Based  on  cross-correlation  peak.  Peak  time  converted  to  range  using  a  typical  freshwater 
sound  speed  of  1430  m/s. 


2.5.2  Current  Divider  Factor 

Using  the  simple  circuit  (Fig.  2),  the  proportionality  factor  a  =  |/rad|/|/|  can  be  com¬ 
puted  using  SPICE.  For  the  resonance  frequencies  212  and  320  Hz,  a  is  0.98  and  0.54, 
respectively. 


2.5.3  Field  Point  Range 

The  m-sequence  used  to  interrogate  the  transducer  is  a  useful  ranging  signal.  By  simply 
cross-correlating  the  drive  voltage  signal  against  the  monitor  hydrophone  signal,  the  direct 
path  arrival  can  be  observed  as  the  first  prominent  peak.  For  a  sample  rate  of  8000  Hz,  the 
sample  period  was  0. 125  ms.  Three  transmissions  from  the  2006  Fake  Washington  exercise 
were  found  to  have  acceptable  drive  current  and  monitor  channel  data.  Table  4  shows  the 
ranges  estimated  for  these  three  transmissions. 


2.5.4  Field  Pressure 

The  monitor  channel  data  autospectra  were  estimated  for  each  corresponding  file  using  the 
Matlab  command 

[G,f]  =  pwelch (x, kaiser (8192) , 0, 8192, 8000, ' onesided' ) ; 
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Figure  9:  An  estimated  monitor  hydrophone  autospectrum,  2006  Lake  Washington  test,  file 
mptest 0  6 . dat. 


This  routine  provides  an  autospectral  estimate  in  units  of  quanta^/Hz.  An  example  au¬ 
tospectrum  is  shown  in  Fig.  9. 

The  ITC  6050C  hydrophone  had  a  terminal  sensitivity  of  —178.2  dB  V^//rPa,  and  the 
A/D  converter  was  set  to  convert  5V  into  2048  quanta,  or  52.2  dB  quanta/V.  The  total 
channel  had  a  transfer  function  of  —126.0  dB  quanta//iPa,  and  was  essentially  frequency 
independent  in  the  circuit  passband.  The  levels  of  the  resonance  peaks  were  extracted  and 
converted  to  in-water  face-of-phone  units  (Table  5). 


pressure 

current 

file 

^,^(212) 

^^^(212) 

Ggg{212) 

G'n(212) 

mptest 04 . dat 

12.8 

138.7 

20.1 

-24.2 

mptest 05 . dat 

12.9 

138.9 

20.2 

-24.1 

mptest 0  6 . dat 

12.6 

138.6 

20.3 

-24.0 

average 

138.7 

-24.1 

amplitude 

8.64  xlOVPa 

6.24  X  10-^  A 

file 

^,^(320) 

^^^(320) 

Ggg{320) 

G'n(320) 

mptest 04 . dat 

-2.1 

123.9 

0.2 

-44.0 

mptest 05 . dat 

-2.0 

124.0 

0.3 

-44.0 

mptest 0  6 . dat 

-2.0 

124.0 

0.3 

-44.0 

average 

123.4 

-44.0 

amplitude 

1.58  xlOVPa 

6.3  X  10"^  A 

Table  5:  Radiated  pressure  and  drive  current  magnitudes  for  the  three  transmissions,  reso¬ 
nances  of  212  Hz  and  320  Hz,  Lake  Washington  2006  test.  The  averages  were  computed  in 
decibel  space.  LFnits  are  dB  re  1  quanta^/Hz  for  Gqq{f),  dB  re  1  /iPa^/Hz  for  Gpp{f),  and 
dBre  1  A^/Hzfor  Gaif). 
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Figure  10:  An  estimated  drive  current  autospectrum,  2006  Lake  Washington  test,  file 
mptest 0  6 . dat. 


resonance  [Hz] 

^2 

5 

5.70  X  10^ 
7.56  X  10^ 

7.55  X  10^ 
8.69  X  10^ 

“average” 

8.12  X  10^ 

Table  6:  Transformation  ratios  calculated  from  the  2006  Lake  Washington  data.  Density 
was  1000  kg/m^,  and  sound  speed  was  1430  m/s.  All  other  values  given  in  the  text. 

2.5.5  Drive  Current 

The  drive  current  data  autospectra  were  similarly  estimated  for  each  corresponding  file 
using  the  Matlab  command 

[G,f]  =  pwelch (x, kaiser (8192) , 0, 8192, 8000, ' onesided' ) ; 

This  routine  provides  an  autospectral  estimate  in  correct  units  of  quanta^/Hz.  An  example 
autospectrum  is  shown  in  Fig.  10. 

The  current  monitoring  circuit  had  a  calibration  of  0.0061  quanta/A.  The  levels  of  the 
autospectra  at  the  resonance  frequencies  are  also  provided  in  Table  5. 


2.5.6  Incorporating  the  Transformation  Ratio 

Using  the  values  of  field  pressure  level,  slant  range,  drive  current  level,  the  local  sound 
speed,  and  coarse  estimates  of  density  and  simple  source  radii,  values  for  the  transformation 
ratio  can  be  calculated  at  resonance  frequencies  of  212  and  320  Hz  (Table  6). 

Using  the  electro-acoustic  transformer,  it  is  now  possible  to  “reflect”  the  two  resonant 
branches  from  the  electrical  equivalent  circuit  to  the  acoustic  stage.  The  “average”  trans- 
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element 

initial  value 

transformed  value 

Clo 

Llo 

Rlo 

Chi 

Lri 

Rhi 

1.08  X  10"^ 
5.22 

111 

2.74  X  10“® 
9.04 

698 

111 

7.92  X  10-1° 

1.93  X  10-® 

181 

1.37  X  10-° 

1.06  X  10-^ 

Table  7:  Lumped  acoustical  equivalent  circuit  values  for  the  resonant  branch  components. 


Figure  1 1 :  SPICE  representation  of  an  ideal  transformer  with  transformation  ratio  S. 


formation  ratio  was  used  in  this  calculation.  This  reflection  is  not  necessary,  but  it  makes 
clear  the  separate  functionality  of  the  acoustic  components.  The  resulting  equivalent  circuit 
is  shown  in  Fig.  7  and  the  transformed  circuit  quantities  in  Table  7. 

SPICE  provides  a  transformer  element,  but  it  is  too  realistic  to  use  here.  The  proper 
“ideal  transformer”  is  shown  in  Fig.  11.  The  “secondary”  consists  of  a  voltage  source 
el  equal  to  S  times  the  primary  voltage,  and  a  “voltage  sensor”  va.  The  voltage  sensor 
supplies  no  additional  voltage,  but  is  a  mechanism  for  measuring  the  current  out  of  the 
secondary.  The  “primary”  consists  of  a  current  source  /I  that  drives  a  current  S  times 
that  flowing  out  of  the  secondary.  The  formal  SPICE  circuit  corresponding  to  the  modified 
equivalent  circuit  Fig.  7  is  shown  in  Fig.  12. 
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3  Tuning 

It  is  advantageous  to  add  a  tuner/transformer  to  a  high-power  piezoceramic  transducer  for 
two  reasons: 

1 .  The  inductance  of  the  tuner/transformer  can  be  used  to  match  the  capacitance  of  the 
piezoceramic  itself,  resulting  in  a  better  power  factor  at  the  power  amplifier,  which 
results  in  more  efficient  transfer  of  power  to  the  transducer; 

2.  The  voltage  on  the  suspension  cable  can  be  kept  within  rated  levels  by  utilizing  a 
tuner/transformer  at  the  underwater  package  that  includes  a  step-up  ratio. 


An  autotransformer  is  generally  used  in  this  application,  because  electrical  isolation  is 
not  required,  and  the  autotransformer  does  not  require  as  massive  a  core  as  a  conventional 
transformer. 


3.1  Single  Resonance  Theory 

Consider  an  example  circuit  (Fig.  13).  At  the  resonance  frequency,  the  untuned  circuit  has 
a  power  factor  of  0.55;  the  admittance  loop  is  shown  in  Fig.  14.  An  inductance  value  of 
50  mH  brings  the  power  factor  to  0.996:  the  resulting  admittance  loop  is  also  shown  in 
Fig.  14. 

Using  the  equation 


P{f)  =  5  v(/)'/(/)  =  w{f)  +  Mf)  (18) 

for  the  average  complex  power,  a  SPICE  calculation  at  the  resonance  frequency  shows  that 
the  ratio  of  real  power  delivered  by  the  source  {W)  to  total  power  \P{f)  \  goes  from  55%  to 
100%,  corresponding  to  the  power  factor.  Delivery  of  power  at  this  frequency  is  optimized 
by  this  inductor. 


3.2  Multiple  Resonance  Theory  and  Wideband  Signals 

The  conventional  approach  presented  in  the  previous  section  is  defined  for  a  system  with 
a  single  resonance  at  a  single  frequency.  It  is  not  obvious,  however,  how  to  extend  those 
results  to  a  system  with  multiple  resonances  and  wideband  signals. 
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Drs.  Birdsall  and  Metzger  (University  of  Michigan)  compute  power  for  the  HX554 
transducers  fabricated  under  the  ATOC  program  as  follows.  Given  time-domain  voltage 
v{t)  and  current  i{t)  at  the  driving  point,  the  power  is 

W  =  ^  v{t)i{t)  dt  (19) 

Based  on  the  Schwarz  inequality,  the  quantity 


will  always  be  less  than  or  equal  to  unity.  The  value  of  unity  is  obtained  if  v(t)  is  linearly 
related  to  i{t).  Furthermore,  for  v{t)  =  Vocos{ut)  and  i{t)  =  Iocos{ut  +  4>),  this  ex¬ 
pression  yields  the  single  frequency  power  factor  of  cos  0.  Eq.  20  is  therefore  a  candidate 
“broadband  power  factor.” 

Consider  a  circuit  (Fig.  15)  involving  two  parallel  branches  to  represent  the  two  reso¬ 
nances  of  the  transducer.  It  includes  a  tuning  inductor  Lp.  Several  values  for  this  tuning 
inductor  are  used:  220  mH,  300  mH,  and  450  mH  (see  section  3.3).  The  additional  values 
provide  coverage  around  the  broadband  power  factor  peak.  The  time  domain  drive  sig¬ 
nal  was  an  m-sequence  with  1023  bits,  a  Q  of  2,  and  carrier  frequencies  of  250,  275,  and 
300  Hz.  These  options  span  the  likely  values  envisioned  for  the  multiport  system.  The  drive 
current  was  computed  by  inverse  Fourier  transformation  of  the  input  admittance  product  (as 
computed  by  SPICE)  and  the  drive  signal  transform  (Table  8). 


Figure  13:  Tuning  example  circuit.  The  shunt  inductor  Lt  is  added  to  compensate  for  the 
reactance  in  the  parallel  combination  of  the  Co  capacitor  and  the  series  R-L-C  branch.  The 
resistor  Rt  is  included  to  model  the  small  resistance  in  the  inductor  windings.  Co  =  10/iF, 
Cl  =  0.1/iF,  Ri  =  100  U,  and  Li  =  5  H.  The  tuning  inductor  is  Lt  =  50  mH,  and 
Rt  =  0.1  U.  The  resonance  frequency  is  ~  225  Hz. 
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Figure  14:  Admittance  loops  for  the  single  resonance  circuit,  with  and  without  the  tuning 
inductor. 


O2 


Figure  15:  Simple  electrical  equivalent  SPICE  model  of  transducer  with  a  parallel  tuning 
inductor.  A  resistor  Rp  is  included  with  the  inductor  to  prevent  a  short  circuit  shunt  path  at 
DC  in  SPICE  analysis.  Rp  is  1  D. 
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carrier  frequency 

Lp 

250  Hz 

275  Hz 

300  Hz 

220  mH 

09.0% 

08.7% 

08.2% 

300  mH 

11.6% 

11.1% 

10.4% 

450  mH 

15.0% 

14.0% 

12.6% 

600  mH 

16.5% 

15.1% 

13.2% 

850  mH 

17.3% 

15.4% 

13.2% 

1000  mH 

17.3% 

15.4% 

13.1% 

1250  mH 

17.2% 

15.1% 

12.8% 

Table  8:  “Broadband  power  factor”  for  the  circuit  of  Fig.  15  and  an  m-sequence  drive 
signal. 


Optimal  tuning  does  not  bring  this  measure  near  100%,  as  was  the  case  in  the  single  fre¬ 
quency  example.  Indeed,  over  the  values  of  tuning  inductance  used,  there  is  little  significant 
change  in  this  measure.  If  the  broadband  power  factor  were  interpreted  as  cos  0,  then  the 
phase  angle  is  generally  85°  to  at  best  80°.  Furthermore,  there  is  a  slight  dependence  on  the 
drive  signal  itself,  a  feature  that  complicates  the  design  of  the  optimal  tuning  inductance. 
Alternatively,  the  measure  is  not  strongly  sensitive  to  these  drive  signals,  so  an  approxi¬ 
mate,  albeit  suboptimal,  tuning  inductance  value  could  be  determined  using  essentially  any 
signal  from  a  given  class  of  target  signal  types. 

This  multiple  resonance  example,  unlike  the  single  one,  shows  that  it  is  neither  simple 
nor  straightforward  to  optimally  tune  a  system  for  a  broadband  signal,  especially  a  system 
with  multiple  resonances.  (On  a  similar  note,  John  Watson  of  APL-UW  [9]  advised  that  a 
computer  analysis  be  performed  over  the  full  band  to  gauge  the  efficacy  of  a  tuner.) 

The  electrical  power  calculated  by  Eq.  19  is  all  real,  i.e.,  “active”  power,  and  does  not 
include  reactive  power. 


3.3  The  1:2  l\iner 

Dr.  Jack  Butler  and  Victoria  Curtis  at  ImageAcoustics,  Inc.,  considered  the  equivalent 
circuit  of  Fig.  15  and  suggested  optimal  values  for  the  parallel  inductor  Lp  based  on  the 
following  argument.  The  curves  of  conductance  and  susceptance  against  frequency  for  the 
untuned  circuit  (i.e.,  without  a  parallel  inductor)  are  shown  in  Fig.  16.  Using  a  design  goal 
of  achieving  the  greatest  number  of  zero-crossings  (or  zero-touchings)  of  the  susceptance 
curve,  they  chose  a  parallel  inductance  value  of  0.296  H.  This  moves  the  susceptance  curve 
down  until  four  points  either  touch  or  cross  the  axis:  see  Fig.  17. 
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Figure  16:  Conductance  and  susceptance  of  the  circuit,  no  tuning  inductor. 


Figure  17:  Conductance  and  susceptance  of  the  circuit,  tuning  shunt  inductor  of  0.296H, 
shunt  resistance  of  1  ohm. 


As  alternatives,  two  other  shunt  inductances  were  considered,  designed  to  cancel  the 
susceptance  at  either  of  the  two  resonances,  in  accord  with  the  simpler  single-resonance 
problem. 

The  first  autotransformer/tuner  was  therefore  specified  to  have  shunt  inductances  of 
220  mH,  300  mH,  and  450  mH.  The  device  was  also  specified  to  have  a  step-up  ratio  of 
primary  to  secondary  voltages  of  1:2. 

The  contract  for  this  device  was  awarded  to  Coiltron,  Inc.  [Tigard,  OR].  Pictures  of  the 
device  when  delivered  are  shown  in  Figs.  18  and  19  .  Concern  regarding  heat  build-up  in 
the  device  led  APL-UW  to  provide  two  thermistors  to  Coiltron  during  the  winding  process 
to  be  inserted  deep  inside  the  windings,  near  the  core,  for  monitoring  during  operation. 
Leads  from  these  thermistors  can  be  seen  in  the  pictures. 

The  housing  for  the  tuner  was  designed  and  fabricated  at  APL-UW  (Figs.  20  and  21). 
The  housing  was  filled  with  oil  to  act  as  a  convective  heat-transfer  coolant  during  operation, 
and  also  to  compensate  hydrostatic  pressure.  Concern  regarding  thermal  expansion  of  the 
oil  when  heated  led  to  the  addition  of  an  oil  exchange  system  (Fig.  22),  which  has  tubing 
to  connect  oil  in  the  housing  with  two  expansion  tanks. 
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Figure  18:  Corner  view  of  the  first  tuner. 


Figure  19:  End  view  of  the  first  tuner. 
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Figure  20:  Housing  for  the  tuner. 


Figure  21:  Top  view  looking  inside  the  tuner  housing,  with  tuner  in  plaee. 
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Figure  22:  Tuner  and  housing  completely  assembled.  Note  the  oil  expansion  system  tubing 
and  expansion  tanks.  Also  visible  on  top  (protected  with  a  bracket)  is  an  over-pressure 
relief  valve. 

3.4  The  1:3  l\iner 


The  1:2  tuner  did  not  provide  enough  output  voltage  (see  section  7.3),  so  a  1:3  tuner  was 
designed.  In  the  SPICE  circuit  for  this  tuner  (Fig.  23),  the  autotransformer  is  modeled  as 
two  magnetically  coupled  inductors,  with  values  Lpi  and  Lp2,  respectively.  There  is  an 
additional  resistance  associated  with  each  inductor  for  modeling  the  resistive  loss  through 
the  windings.  The  magnetic  coupling  factor  used  in  the  calculations  was  0.98. 


The  equations  governing  the  autotransformer  are  as  follows.  The  primary  voltage  vi 
and  primary  current  ii  are  related  to  the  secondary  voltage  V2  and  secondary  current  *2  by 


vi  =  jujLpi{ii  +  i2)  +  jujMi2  (21) 

V2  =  jujLp2i2  +  jujLpi{ii  +  is)  +  jo;  Mil  (22) 


where  M  is  the  mutual  inductance,  defined  according  to  the  magnetic  coupling  factor  K  as 


K 


M 


^  Lp\Lp2 


(23) 


The  simplest  analysis  assumes  perfect  coupling  (i.e.,iF  =  1,  and  therefore  M  =  ^LpiLp2) 
and  an  open  circuit  secondary.  In  this  case,  is  =  0  and  the  equations  above  become 
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Figure  23:  Simple  electrical  equivalent  SPICE  circuit  involving  the  autotransformer. 
Solving  for  the  output/input  ratio,  one  obtains 


Vl 


(26) 


This  approximation  is  valid  when  the  inductive  reactance  is  much  greater  than  the  series 
resistance  in  each  winding,  and  the  secondary  load  is  large  compared  to  the  inductive  reac¬ 
tance  seen  across  the  entire  secondary.  This  will  in  general  not  be  the  case  when  operating 
near  the  tuned  frequency,  but  this  approximation  does  provide  a  useful  crude  formula  to 
verify  the  step-up  ratio  of  the  autotransformer. 


(The  choice  of  lower  and  upper  inductance  values  —  here  Lpi  and  Lp2  —  to  meet  a 
required  step-up  ratio  and  overall  inductance  involves  more  complexity  than  is  represented 
in  Eq.  26  [pers.  comm.,  Richard  Gettmann,  13  July  2009].  Therefore,  Eq.  26  should  be 
interpreted  as  a  rough  guide  for  choosing  the  two  inductance  values.  It  remains  true  that 
Eq.  26  provides  a  good  model  for  the  output  of  the  corresponding  SPICE  circuit. 


The  design  goal  was  to  achieve  a  step-up  from  the  primary  to  the  secondary  of  about  a 
factor  of  3,  and  to  adjust  the  susceptance  curve  about  zero.  This  was  done  by  trial- and-error 
at  each  of  the  two  resonance  frequencies,  thereby  determining  the  “optimal”  solution  for 
each.  The  middle  option  was  devised  as  was  done  for  the  1:2  tuner,  by  determining  those 
values  that  caused  four  crossings  or  touchings  of  the  susceptance  curve  with  the  zero  axis. 
The  resulting  parameter  values  are  shown  in  Table  9. 

The  admittance  loops  for  these  three  cases  are  shown  in  Eigs.  24,  25,  and  26,  and  it 
can  be  seen  for  the  lower  and  upper  resonances  that  the  proper  choice  of  autotransformer 
inductors  centers  the  corresponding  admittance  loop  on  the  zero  axis.  The  “middle”  tuning 
configuration  centers  neither  on  the  zero  axis. 
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Lp\ 

Lp2 

low  freq 
mid  freq 
hi  freq 

50  mH 
33  mH 
22  mH 

210  mH 
130  mH 
90  mH 

Table  9:  Upper  and  lower  inductor  values  for  tuning  the  circuit,  at  three  different  frequen¬ 
cies,  1:3  step-up. 

The  corresponding  conductance  and  susceptance  curves  for  these  three  cases  are  shown 
in  Figs.  27,  28,  and  29.  The  “middle”  tuning  can  be  verified  here  as  that  configuration  that 
causes  four  crossings  or  touchings. 

As  an  additional  verification,  the  tuner  output  (autotransformer  secondary)  is  plotted 
versus  frequency  for  all  three  tuning  configurations  in  Figs.  30,  31,  and  32.  The  input 
voltage  in  all  cases  was  uniform  across  all  frequencies  at  1  -f  Oj.  These  figures  indicate  that 
the  autotransformer  secondary  output  is  load  dependent  for  these  load  magnitudes,  but  that 
the  output-input  ratio  is  close  to  the  design  value  of  3  for  all  frequencies  considered. 
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Figure  24:  Admittance  loop,  tuned  for  the  lower  (212  Hz)  resonance  frequency. 


Figure  25:  Admittance  loop,  tuned  between  the  resonances. 


Figure  26:  Admittance  loop,  tuned  for  the  upper  (320  Hz)  resonance  frequency. 
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Figure  27:  Conductance  and  susceptance  curves,  tuned  for  the  lower  (212  Hz)  resonance 
frequency. 


Figure  28:  Conductance  and  susceptance  curves,  tuned  between  the  resonances. 


Figure  29:  Conductance  and  susceptance  curves,  tuned  for  the  upper  (320  Hz)  resonance 
frequency. 
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Figure  30:  Tuner  voltage  across  the  secondary.  The  input  voltage  was  unity.  Circuit  tuned 
for  the  lower  (212  Hz)  resonance  frequency. 


Figure  31:  Voltage  across  the  secondary.  The  input  voltage  was  unity.  Circuit  tuned  be¬ 
tween  the  resonances. 


Figure  32:  Voltage  across  the  secondary.  The  input  voltage  was  unity.  Circuit  tuned  for  the 
upper  (320  Hz)  resonance  frequency. 
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4  The  Multiport  System 


An  essential  system  bloek  diagram  of  the  APL-UW  system  (Fig.  33)  shows  the  underwater 
subsystem  eomponents: 

•  the  MP200/TR1446  with  the  tuner 

•  a  monitor  hydrophone 

•  a  eustom  telemetry  bottle  with  a  fibre-optie  interfaee 

•  a  depth  sensor 

•  a  SeaBattery 


Figure  33:  Block  diagram  of  essential  components  for  the  complete  multiport  system. 
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The  ship-side  equipment  consists  of  a  signal  delivery  rack,  which  includes  a  80686 
PC  running  Fedora  Core  8  Linux,  a  Spectrum  Instruments  TM-4  GPS  receiver,  a  custom 
interface  unit,  and  a  Krohn-Hite  filter;  a  power  amplifier;  a  custom  winch;  and  a  tracking 
system.  The  original  power  amplifier  was  a  Ling.  An  L50  amplifier  was  procured  from 
Instruments,  Inc.,  for  the  Philippine  Sea  experiments. 

The  underwater  package  is  suspended  from  an  “electro-opto-mechanical”  cable.  Teleme¬ 
try  information  utilizes  the  fibre,  and  transducer  power  utilizes  the  copper  conductor.  In 
practice,  the  underwater  package  is  deployed  and  recovered  by  the  custom  winch. 


An  ITC  6050C  hydrophone  will  normally  be 
deployed  on  the  transducer  axis  at  a  known,  fixed 
distance  (typically  20  m)  above  the  acoustic  cen¬ 
ter  of  the  main  tube.  The  raw  time-domain  acous¬ 
tic  signal  is  digitized  in  the  telemetry  bottle,  mul¬ 
tiplexed  onto  the  optical  fibre,  demultiplexed  and 
reconstructed  back  to  an  analog  signal  within  the 
shipboard  fibre  interface,  and  supplied  to  the  sig¬ 
nal  computer  for  recording.  Initial  measurements 
of  this  channel  (not  including  the  hydrophone)  in¬ 
dicate  that  the  channel  has  a  passband  level  of  —5.5 
dB  between  about  7  Hz  and  5000  Hz. 

Under  the  assumption  that  the  acoustical  per¬ 
formance  of  the  MP200/TR1446  itself  would  be 
compromised  if  anything  was  attached  to  it,  the 
mechanical  configuration  hung  the  auxiliary  com¬ 
ponents  (SeaBattery,  telemetry  bottle,  and  tuner) 
above  the  MP200/TR1446  (Fig.  34).  In  this  de¬ 
sign,  the  MP200/TRI446  itself  would  be  stowed 
on-board  deck  horizontally. 


Subsequently,  in  consultation  with  ImageA-  piguj-e  34:  Schematic  of  the  initial 
coustics,  it  was  determined  that  the  configuration  mechanical  configuration, 
in  Fig.  34  would  mass-load  the  fluid  in  the  main 
tube  and  alter  the  resonant  and  radiation  properties 
of  the  device.  (It  was  also  realized  that  the  chain  of 
components  may  be  extremely  difficult  to  deploy.) 


The  mechanical  design  evolved  so  that  the  auxiliary  components  were  relocated  to  the 
equator  of  the  MP200/TR1446  and  placed  on  opposite  sides  for  balance  (Fig.  35).  An 
auxiliary  structure  was  designed  that  both  supported  the  auxiliary  components  and  also 
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supported  the  MP200/TR1446  in  an  upright  position  while  on  deck.  Although  massive,  the 
SeaBattery  and  the  tuner  were  small  relative  to  a  wavelength,  and  therefore  their  modifi¬ 
cation  to  the  radiated  acoustical  field  should  be  small  as  long  as  they  are  away  from  the 
tube  openings.  The  support  structure  itself  was  fabricated  from  free-flooded  steel  tubing, 
which  also  should  have  a  negligible  acoustical  cross-section.  No  modeling  has  been  done 
to  quantify  the  effect  of  these  auxiliary  components  on  the  radiated  acoustical  field. 

With  the  design  complete,  it  was  realized  that  the  entire  system  would  be  too  tall  to  be 
lifted  from  the  deck  using  a  chain  bridle  attached  to  the  bolt  eyes  fixed  into  the  top  of  the 
main  tube.  A  tentative  plan  was  devised  to  lift  from  a  hard  point  welded  onto  the  main  tube 
exactly  at  the  top  of  the  secondary  tube,  i.e.,  lifting  the  entire  system  from  one  “shoulder.” 
This  plan  introduced  a  further  problem:  the  system  would  hang  at  an  angle  significantly  off 
vertical.  Further  calculations  showed,  however,  that  a  lifting  point  located  at  the  center  of 
a  cross-bar  bolted  into  the  mouth  of  the  main  tube  several  centimeters  down  from  the  top 
would  provide  enough  clearance.  The  cross-bar  width  was  much  less  than  a  wavelength 
and  overall  cross-sectional  area  much  less  than  the  area  of  the  main  tube  opening,  so  this 
bar  is  not  expected  to  modify  the  resonant  or  radiation  patterns  significantly.  In  the  final 
analysis,  the  system  was  calculated  to  hang  about  2.6°  off  vertical. 
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Figure  35:  Engineering  drawings,  MP200/TR1446  assembly. 
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Figure  36:  System  fully  assembled  on  the  deck  of  the  M/V  SeaHorse,  2009  Lake  Washing¬ 
ton  test. 
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Figure  37:  Deck  preparations,  2009  Lake  Washington  test. 
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Figure  38:  Deployment  from  the  M/V  SeaHorse,  January  2009  Lake  Washington  test. 
System  fully  assembled  ineluding  sea  battery,  telemetry  bottle,  and  tuner.  Deployment  and 
reeovery  utilized  the  barge’s  main  erane.  The  lifting  strap  is  attaehed  to  a  hard  point  on  the 
transdueer;  this  was  added  based  on  an  early  lifting  design.  Onee  in  the  water,  the  system 
is  transferred  to  the  “eleetro-opto-meehanieal”  suspension  eable,  and  is  raised  and  lowered 
in  the  water  via  the  APL-UW  wineh  and  suspension  eable.  In  this  pieture,  the  suspension 
eable  is  only  under  mild  tension.  The  bloek  is  also  suspended  from  the  main  erane.  The 
eable  married  for  a  short  distanee  up  the  suspension  eable  is  the  monitor  hydrophone  eable 
(hydrophone  not  seen). 
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5  System  Model 


5.1  Introduction 

The  simple  equivalent  circuit  model  (section  2),  which  models  only  the  electrical  character¬ 
istics  of  the  transducer,  can  be  expanded  to  include  all  system  components,  from  the  input 
data  file  to  the  output  radiated  far  field  pressure.  This  provides  a  complete  end-to-end  mod¬ 
eling  capacity.  A  block  diagram  of  the  system  is  shown  in  Fig.  39.  Such  a  model  requires 
the  additional  “transformation  ratio”  in  the  equivalent  circuit  model  to  transform  from  the 
electrical  stage  to  the  acoustical  lumped  equivalent  circuit,  in  addition  to  elements  for  all 
other  circuit  components,  such  as  the  computer  D/A,  the  power  amplifier,  etc. 


5.1.1  D/A  Converter 


The  D/A  converter  is  one  subsystem  on  the  National  Instruments  NI  PCI  607  IE  DAQ 
board.  This  is  a  12-bit  converter.  The  board  uses  the  COMEDI  driver,  version  7.76.  The 
conversion  from  input  quanta  to  output  voltage  is  set  to  10  V/2048  quanta,  or  4.88  mV  per 
quanta. 

Modeling  this  element  requires  only  a  simple  scaling  operation,  where  the  scaling  factor 
is  4.88  mV/quanta.  This  can  be  achieved  very  simply  in  a  SPICE  circuit  with  a  2-port 
network  that  drives  a  scaled  version  of  the  input  “voltage”  to  the  output  port  (Eig.  40).  In 
this  case,  the  input  signal  is  not  measured  in  volts  but  in  D/A  quanta. 


Eigure  39:  Block  diagram  of  the  entire  system.  The  individual  components  are  modeled  as 
SPICE  “subcircuits.” 
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v2  =  A  vl 


Figure  40:  SPICE  model  for  an  ideal  D/A  converter.  For  this  system,  scale  factor  A  is 
4.88  mV/quanta.  This  is  a  generic  model  for  all  simple  signal  “scaling”  operations  that  do 
not  require  fidelity  of  the  input  and  output  currents  and  impedances. 


drive  voltage 
[V  RMS] 

output  voltage 
[V  RMS] 

ratio 

621 

258.8 

2.74 

708 

258.4 

799 

258.6 

3.44 

889 

258.4 

Table  10:  Input  and  output  measurements  for  a  300-Hz  sine  wave,  Instruments,  Inc.,  L50 
power  amplifier.  Source  setting  is  voltage  source  with  tap  1200/5.4. 


5.1.2  Power  Amplifier 

For  simplicity,  the  power  amplifier  is  modeled  as  a  simple  scaling  factor,  and  therefore  uses 
the  same  circuit  as  shown  in  Fig.  40. 

The  target  power  amplifier  for  at-sea  experiments  was  an  Instruments,  Inc.,  [San  Diego, 
CA]  F50,  modified  for  operation  at  low  frequencies.  The  transfer  function  of  the  amplifier 
was  measured  at  APF-UW  using  a  sine  wave  generator  for  input,  and  a  500f2  resistor  for 
a  dummy  load.  Input  signal  RMS  level  was  measured  on  a  multimeter:  amplifier  output 
voltage  RMS  level  was  measured  by  a  custom  monitoring  circuit.  To  a  good  approximation, 
the  amplifier  acts  as  a  signal  scaling  element  with  a  scale  factor  of  258  (Table  10). 


5.1.3  Low-Pass  Filter 


To  avoid  a  175-Hz  4-pole  low-pass  Bessel  filter  built  into  the  custom  interface  unit,  the  raw 
D/A  output  from  the  interface  unit  was  used.  This  raw  signal  requires  low-pass  filtering 
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Figure  41:  Circuit  model  for  the  low -pass  filter. 


element 

value 

comment 

1  n 

small,  to  minimize  network  loss 

C'l 

31.8 /iF 

to  achieve  a  comer  at  5000  Hz 

Table  11:  Circuit  parameters  for  the  low  pass  filter  model. 


to  smooth  out  the  staircase  waveform  output  from  the  D/A.  This  was  accomplished  with 
a  Krohn-Hite  3343  rack-mount  filter.  (A  typical  low-pass  filter  corner  setting  is  600  Hz 
on  the  Krohn-Hite  filter.)  A  simple  1-pole  RC  (resistor-capacitor)  network  was  used  to 
crudely  model  the  low-pass  behavior  of  the  filter  (Fig.  41).  Circuit  element  values  are 
provided  in  Table  11. 


5.1.4  Cable 

Fred  Karig  (APL-UW)  provided  the  cable  parameters  (Table  12).  Fig.  42  shows  the  elec¬ 
trical  model  for  the  suspension  cable. 


5.2  Acoustic  Radiation:  the  “Double  Doublet”  Model 

A  simple  model  for  the  acoustic  radiation  from  a  single  tube  (section  2.5)  models  the  tube 
as  an  acoustic  doublet.  Extending  this  concept  to  the  MP200/TR1446  produces  a  “double 
doublet”  (Fig.  43).  The  radiation  from  each  tube  opening  —  the  piston  in  the  end  of  an 
unfianged  (or  partially  flanged)  pipe  —  is  idealized  as  a  monopole  in  free  space.  Let  the 


element 

value 

comment 

Rcable 

3.15  D 

1.5  D/IOOO  ft  from  the  manufacturer 

C cable 

10.5  /xF 

25  pF/ft ,  using  RG-8  specification 

Table  12:  Circuit  parameters  for  the  cable.  Cable  length  is  4200  ft. 
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Figure  42:  Circuit  model  for  the  cable. 


acoustic  particle  velocity  of  the  longer  tube  (lower  frequency)  be  Mlo  and  that  from  the 
shorter  tube  (higher  frequency)  be  Mhi-  The  Euclidean  vector  from  the  monopole  to  a  field 
point  is  rLo/Hi+ /- ,  where  the  subscripted  sign  refers  to  the  vertical  position  of  the  monopole. 

The  total  far  field  pressure  then  utilizes  the  free  space  Green  function  from  each  mono¬ 
pole,  and  for  a  single  frequency  is  given  by 


p(r)  =  n 


^ikrLO+  pikrHi+ 


Jkrni- 


LO 


47rr 


n 


HI 


LO+ 


Anr 


n 


HI 


HI+ 


Atithi- 


n 


LO 


^ikrto- 


Atitlo- 


(27) 


where  the  source  level  factors  are  given  by  (see  Eq.  11) 


^LO,HI  —  '^0-LO,HlPC^'^LO, 


HI 


(28) 


5.3  Model  Implementation 

The  system  model  and  the  adjunct  radiation  routines  were  combined  into  a  C-f-f  program. 
The  code  generated  the  SPICE  netlisting  for  the  model,  and  initiated  a  call  to  SPICE  to 
evaluate  the  circuit  at  multiple  frequencies.  A  typical  netlist  is  shown  in  Appendix  B.  The 
code  then  parsed  the  resulting  SPICE  output  file  to  retrieve  the  transducer  fluid  velocities 
Ulo  and  uhi-  The  code  then  used  Eqs.  27  and  28  to  compute  the  complex  acoustic  pres¬ 
sure  response  at  each  frequency  at  the  desired  field  point.  The  code  next  multiplied  the 
input  signal  spectrum  by  the  complex  acoustic  pressure  response  at  each  frequency,  and 
inverse  Eourier  transformed  the  result  to  a  real  time-domain  signal.  User-adjustable  vari¬ 
ables  include  the  sound  speed,  the  field  point  {x,  y,  z)  in  a  coordinate  system  centered  at 
the  transducer,  and  the  file  name  of  the  desired  drive  signal  in  .  wav  format. 

A  variety  of  predictions  can  be  made  with  this  model  and  code.  Some  examples  follow. 
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5.4  Performance  Validation 

Because  the  RMS  acoustic  pressure  at  [30.8,  0,  0]  was  used  to  relate  ulo,hi  to  the  transfor¬ 
mation  ratio  S  via  Eqs.  17  and  8  via  Eqs.  9  and  1 1  and  a  simplified  version  of  the  “double 
doublet”  propagation  model,  the  reverse  calculation  must  hold:  a  sinusoidal  input  at  212  Hz 
should  yield  the  measured  RMS  value  at  the  field  point.  The  measured  RMS  pressure  was 
8.64  X  10®  /iPa  at  RMS  input  current  of  6.24  x  10^  A  (Table  5).  (This  current  was  measured 
at  the  top  of  the  cable,  but  this  will  not  introduce  a  significant  error.) 

The  SPICE  netlisting  was  altered  to  introduce  a  current  sensing  voltage  source  between 
the  tuner  and  the  transducer.  Simple  SPICE  runs  at  the  single  frequency  212  Hz  showed 
that  a  power  amplifier  factor  of  8.83  produced  the  target  RMS  current  for  a  Vm  input  with 
amplitude  100. 

Next,  the  simulator  code  was  modified  to  force  the  power  amplifier  scale  factor  to  be 
8.83,  and  the  code  exercised  on  an  input  file  containing  a  sinusoid  of  frequency  212  Hz  and 
amplitude  100  quanta.  The  output  file  at  the  field  point  [30.8,  0,  0]  was  a  simple  sinusoid, 
and  was  read  into  Matlab  where  it  was  found  to  have  an  RMS  value  of  7.09  x  10®,  which  is 
in  error  of  1.7  dB.  This  is  adequate  verification,  considering  all  the  approximations  made 
in  this  simple  transducer  model. 


5.4.1  Broadband  Beampattern 

The  MP200/TR1446  is  known  to  have  a  beampattern  with  non-negligible  deviations  from 
an  omnidirectional  pattern,  particularly  at  the  upper  resonance  frequency.  An  interesting 
relevant  broadband  computation  is  therefore  the  “beampattern”  for  an  input  m-sequence. 
Input  signal  file  m2  84  .  LW .  eq .  wav  was  used  for  this  calculation.  This  file  contains  an 
m-sequence  with  law  2033,  Q  of  2,  carrier  of  284.166666. . .  Hz,  and  a  sample  rate  of 
8525  Hz.  The  amplitude  was  170  quanta.  The  signal  had  been  pre-equalized  using  the 
predicted  equalizer  filter  (see  section  6).  A  Matlab  script  was  used  to  exercise  the  simulator 
at  angles  from  0°  (i.e.,  broadside)  to  90°  (i.e.,  endfire)  at  a  slant  range  of  100  m.  At  each 
angle,  the  script  read  the  output  time-domain  file  (in  units  of  /iPa)  and  computed  the  RMS 
value.  The  normalized  beampattern  is  defined  here  as 

Results  show  a  drop  of  only  a  4  dB  in  the  broadband  RMS  levels  from  broadside  to  endfire 
(Fig.  44). 
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Figure  43:  Geometry  for  the  acoustic  field  computation.  The  radiated  field  is  assumed  to 
be  axially  symmetric  around  the  transducer  long  axis.  The  axial  field  position  at  (0,  0,  20) 
represents  the  monitor  phone  attached  to  the  suspension  cable. 


angle  off  broadside  [degree] 


Figure  44:  Theoretical  broadband  beampattern,  based  on  a  pre-equalized  Q=2  m-sequence 
with  carrier  at  ~  284  Hz. 
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5.4.2  Transmit  Response 


The  full  system  model  can  be  used  to  compute  a  transmit  response  similar  to  the  standard 
transmit  voltage  response  (TVR)  used  in  transducer  design.  In  this  case,  the  transfer  func¬ 
tion  is  from  the  input  .  wav  file  (in  units  of  “quanta”)  to  output  far  field  radiated  pressure 
corrected  back  (via  simple  spherical  spreading)  to  1  m. 

As  an  example,  Fig.  45  shows  the  calculation  for  the  entire  system  model  for  a  far  field 
broadside  direction.  This  computation  was  performed  using  a  Matlab  script  that  looped 
through  all  frequencies;  at  each  frequency,  the  script  called  a  program  to  make  a  .  wav  file 
with  5  s  of  amplitude  100  sinusoid  at  the  target  frequency,  called  the  simulator  to  propagate 
the  field  to  the  position  [100,0,0],  and  then  read  the  resulting  output  file  (which  was  a 
sinusoid  in  units  of  /rPa)  and  computed  its  range-corrected  RMS  value.  The  response  is 
this  value  normalized  by  the  RMS  amplitude  of  the  input  sinusoid,  which  is  70.71  quanta. 


Figure  45:  Theoretical  transmit  response  based  on  the  system  model.  Field  location  was 
[100,  0,  0],  i.e.,  100  m  at  broadside.  Input  units  are  quanta,  i.e.,  the  exact  .  wav  file  values 
supplied  to  the  D/A  converter. 
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5.4.3  Broadside  to  Monitor  Hydrophone  Scaling 

Transmit  level  monitoring  during  actual  at-sea  deployments  will  have  only  the  monitor 
hydrophone  available.  Thus,  it  is  desirable  to  devise  a  relationship  between  the  broadband 
level  measured  at  the  location  of  the  monitor  hydrophone  ([0,0,21])  and  the  broadside 
broadband  source  level  corrected  to  1  m.  This  is  an  easy  computation  to  perform  with 
the  model.  As  an  example,  using  the  input  file  m2  8  4  .  LW .  eq .  wav  and  a  sound  speed 
of  1480  m/s  yields  an  RMS  pressure  of  9.53  x  10®  /iPa  at  [100, 0, 0]  and  2.89  x  10^  at 
[0,  0,  21].  The  ratio  l^^Pfieid.RMs  to  Pmon.RMs  (range-corrected  field  RMS  value  to  monitor 
RMS  face-of-phone  value)  is  33.32,  or  30.5  dB. 
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6  Signal  Equalization 

6.1  Introduction 

The  far  field  transfer  funetion  of  the  system  is  not  flat  aeross  the  two  resonanees,  and  this 
introduees  eonsiderable  speetral  shaping,  typieally  undesirable,  to  any  transmitted  signal. 
To  ameliorate  this  undesirable  shaping,  an  “equalizer”  was  eonstrueted  by  engineers  at 
Massa  Produets,  Ine.  (Fig.  46).  This  unit  was  loaned  to  APL-UW  in  early  2008  under  the 
terms  of  the  tuner  design  eontraet  with  Image Aeousties,  Ine. 

The  transfer  funetion  of  the  Massa  equalizer  was  measured  on  an  HP3589A  speetrum 
analyzer  (Figs.  47  and  48).  This  measurement  required  the  use  of  an  auxiliary  80486  DOS 
maehine  eonfigured  with  a  National  Instruments  GPIB  board  to  eapture  the  traee  to  a  DOS 
ASCII  file  (Table  13).  The  Massa  equalizer  would  typically  be  inserted  between  the  drive 
signal  computer  and  the  power  amplifier  (Fig.  49). 

The  use  of  the  Massa  equalizer  has  the  advantage  of  simplicity:  compensation  for  the 
spectral  characteristics  of  the  transducer  is  achieved  by  simply  adding  a  hardware  compo- 


Figure  46:  The  Massa  equalizer  disassembled. 
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frequency  [Hz] 


Figure  47:  Magnitude  response  of  Massa  equalizer  as  measured  on  an  HP3589  spectrum 
analyzer.  Resolution  bandwidth  1.1  Hz. 


Figure  48:  Phase  response  of  Massa  equalizer  as  measured  on  an  HP3589  spectrum  ana¬ 
lyzer.  Resolution  bandwidth  1 . 1  Hz. 


Figure  49:  Schematic  block  diagram  involving  the  Massa  equalizer.  (Tuner  not  shown). 
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command 

purpose 

ibic 

ibdev  01901310 

ibwrt  "DISPiRES  Dl" 
ibwrt  "FORMiDATA  ASC,7" 
ibwrt  "TRAC: DATA?" 

ibrdf  "FOO.DAT" 

ibonl  0 

starts  the  interactive  session 
opens  a  device,  connects  to  the  analyzer 
move  trace  in  Dl  register  to  trace  A 
set  data  format  to  ASCII,  7  sig.  digits 
queues  trace  A  data  for  retrieval 
transfers  data  from  analyzer  to  DOS  file 
closes  the  device 

Table  13:  AT-GPIB  commands  for  retrieving  trace  data  from  the  HP3589A. 


nent  into  the  signal  circuit.  The  drawback  of  this  approach  is  that  the  model  for  the  entire 
system  now  requires  a  model  for  the  Massa  equalizer,  which  was  not  provided. 

The  alternative  described  below  is  to  construct  a  “pre-equalizer”  filter  to  pre-filter  the 
drive  signal.  This  filter  is  implemented  in  the  frequency  domain.  The  advantage  is  that 
the  system  model  developed  previously  need  not  be  modified,  and,  further,  that  the  filter 
can  be  made  to  compensate  nearly  exactly  for  the  hardware  response  (when  SNR  is  high.) 
The  disadvantage,  primarily  for  the  frequency  domain  filter,  is  that  streaming  operation  of 
an  original  waveform  through  the  filter  is  not  possible  without  recourse  to  a  more  complex 
signal  processing  algorithm  (such  as  an  overlap-and-save  methodology). 


6.2  Theory 

The  equalizer  filter  can  be  developed  using  concepts  from  Wiener  filter  theory,  where  a  sig¬ 
nal  X,  after  possibly  passing  through  some  signal  conditioning,  is  received  in  the  presence 
of  additive  noise  N  (Fig.  50).  Let  the  acoustic  field  at  the  receiver  be 

P{r)  =  XDGH{r).  (30) 


Therefore,  the  received  signal  is 


Y{r)  =  P{r)K  +  NK.  (31) 

The  optimal  filter  is  chosen  to  minimize  the  expected  error  between  F(r)  and  P(r). 
This  gives 

^  A\D\\HG{r)ySxx 

*  mHGir)\^Sxx  +  SxN  ^  ^ 

where  Sxx  is  the  spectrum  of  the  original  input  signal,  is  the  autospectrum  of  the 
additive  noise,  and  A  is  a  user- adjustable  scaling  factor.  Using  simply  Sxx  — ^  \X\‘^  and 
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noise 


Figure  50:  Signal  flow  diagram  for  a  typical  post-processing  equalizer  filter.  D  represents 
the  D/A  converter,  the  power  amplifier,  and  the  step-up  part  of  the  tuner.  H  represents  the 
transducer,  which  is  shown  with  two  outputs,  one  for  each  pair  of  radiating  tube  openings. 
G  is  the  Green  function  from  both  pairs  of  openings  to  a  field  point.  K  is  a  filter  applied 
to  the  data  received  at  a  sensor  at  the  field  point.  X  is  the  input  signal,  P  is  the  acoustic 
signal  at  the  receiver.  Y  is  the  ultimate  equalized  output  signal. 


Eq.  30,  one  obtains 

P-  ^DP-{v) 

P*  ~  \-r>(  \\2  \  C 

|P(r)|^  -f  Snn 

This  is  the  equalizer  transfer  function. 

In  Wiener  filter  theory,  is  computed  from  actual  measurements  of  P(r)  and  Smn- 
Because  the  filter  may  depend  on  direction,  the  choice  of  r  as  the  far  field  broadside  di¬ 
rection  would  be  appropriate,  but  it  was  not  expected  that  measurements  at  such  a  loca¬ 
tion  would  be  available  during  either  calibration  or  actual  field  exercises.  Therefore,  a 
suboptimal  approach  was  defined  as  follows.  The  quantity  P{r)  was  computed  using  the 
SPICE/C-f-f  code,  and  the  denominator  in  Eq.  33  was  replaced  by 

|P(r)p  +  P^^^|P(r)|2(l  +  e)  (34) 

where  e  is  a  frequency-independent  value  that  determines  the  contribution  of  the  prediction 
P(r)  to  the  filter,  e  is  essentially  the  ratio  of  the  noise  to  the  signal.  Eor  e  much  greater  than 
1 ,  the  output  is  dominated  by  the  additive  noise,  and  the  filter  mitigates  the  noise  in  the 
output  by  strongly  attenuating  components  in  X.  Eor  small  e,  the  output  is  dominated  by 
the  signal,  and  the  filter  essentially  deconvolves  out  the  factor  DHG{r).  The  choice  here 
for  e  is  arbitrary  in  this  suboptimal  case. 

All  the  filters  in  Eig.  50  are  linear  operators,  and  therefore  can  be  applied  in  any  order. 
The  equalizer  filter  is  therefore  moved  from  its  “post-acquisition”  position  in  Eig.  50  to  a 
position  prior  to  the  signal  transmission  equipment  in  Eig.  51.  An  off-line  implementation 
of  this  arrangement  is  shown  in  Eig.  52. 

To  demonstrate  the  performance  of  the  equalizer,  a  filter  was  developed  for  the  system 
model  described  previously:  this  filter  is  known  as  file  eqOl  .nc,  with  amplitude  and 
phase  characteristics  shown  in  Eig.  53,  and  e  =  1000.  Note  the  two  deep  nulls  at  the  two 
transducer  resonance  frequencies,  and  the  associated  phase  corrections. 
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Figure  51:  Similar  to  Fig.  50,  but  with  the  equalizer  filter  placed  before  the  D/A  converter. 


Figure  52:  Off-line  implementation  of  pre-equalization,  based  on  Fig.  51. 


An  m-sequence  with  carrier  frequency  300  Hz  and  a  sample  rate  of  9600  Hz  was  pro¬ 
cessed  through  the  system  simulator  and  projected  to  a  location  of  [100, 0, 0],  i.e.,  a  broad¬ 
side  position  at  range  100  m.  The  signal  at  that  location  was  pulse-compressed  using  a 
standard  algorithm:  the  result  is  shown  in  the  top  panel  of  Fig.  54.  This  procedure  was  re¬ 
peated  for  the  same  m-sequence  pre-filtered  by  the  transfer  function  in  file  eqO  1 .  nc:  the 
result  is  shown  in  the  bottom  panel  of  Fig.  54.  Without  equalization,  the  single  pulse,  which 
would  normally  result  from  pulse  compression,  is  distorted  into  multiple  trailing  pulses,  but 
with  equalization  the  pulse-compressed  output  contains  a  single  prominent  pulse. 

Ideally,  the  best  filter  would  use  the  actual  drive  signal  and  actual  field  location  mea¬ 
surements. 


Figure  53:  Equalizer  transfer  function,  file  eqO  1 .  nc 
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Figure  54:  Pulse-compressed  output  without  and  with  equalization.  Simulated  output  at 
location  [100,0,0]:  direct  path  propagation  delay  at  1480  m/s  is  ^  68  ms. 
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Figure  55:  Pulse  compressed  noise  signatures,  2009  Lake  Washington  test,  from  monitor 
hydrophone.  Gain  was  code  5. 

6.3  In-water  Tests 

The  pulse-compressed  signal  from  the  monitor  hydrophone  acquired  during  a  2009  Lake 
Washington  test  is  shown  in  Fig.  55.  The  record  shows  the  direct  path  arrival,  and  the 
surface  reflected  path  arrival.  The  MP200/TR1446  was  at  a  depth  of  30  m.  A  comparison  of 
the  acquired  autospectrum  and  a  predicted  “signal-only”  autospectrum  is  shown  in  Fig.  56: 
this  computation  includes  the  surface  reflected  arrival.  There  is  reasonable  correspondence 
between  the  predicted  spectral  shape  and  that  visible  above  the  background  ambient  noise, 
which  was  filled  with  harmonics  from  the  ship’s  generators.  Note  that  the  nulls  in  the 
predicted  spectrum  do  not  fall  exactly  on  the  nulls  observed  in  the  radiated  spectrum.  This 
suggests  that  the  equalizer  based  on  the  model  is  not  quite  accurate  enough. 


6.4  Summary 

Pre-equalization  based  solely  on  the  system  model  produced  reasonable  results  during  in¬ 
water  tests.  Sharp  resonant  peaks  are  observed  both  in  the  model  and  in  the  actual  device; 
the  sharpness  of  these  peaks,  however,  makes  accurate  compensation  difficult,  as  even  a 
slight  mis-tuning  will  place  spectral  minima  off  actual  resonance.  Such  mis-tuning  can  be 
seen  in  the  actual  radiated  signals. 

Nevertheless,  for  the  m-sequences  used  here,  the  pre-equalizer  produces  strong,  well¬ 
shaped  pulses  (after  pulse  compression)  with  very  low  sidelobes. 
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Comparison:  Real  and  Simulated  Data  Type=monitor 


Figure  56:  Real  and  simulated  radiated  spectra,  2009  Lake  Washington  test  from  monitor 
hydrophone. 
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Pre-equalizing  a  very  long  original  waveform,  or  a  streaming  signal,  would  require  a 
much  more  advanced  algorithm.  For  such  a  case  a  post-signal  generator  hardware  filter, 
i.e.,  the  Massa  equalizer,  is  probably  simpler. 
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7  The  Philippine  Sea  Engineering  Test  /  Pilot  Study 

7.1  Summary 

For  several  deeades,  the  Office  of  Naval  Research  has  sponsored  an  international  consor¬ 
tium  of  scientists  to  investigate  deterministic  and  stochastic  acoustic  propagation  at  low 
frequencies  over  long  ranges  in  the  deep  ocean.  Most  of  these  “blue  water”  experiments 
have  been  conducted  in  the  central  North  Pacific,  where  oceanographic  processes  are  rel¬ 
atively  benign.  The  natural  progression  for  these  studies  is  to  determine  whether  and  to 
what  extent  the  models  and  predictions  developed  during  these  efforts  apply  in  a  region 
with  more  vigorous  oceanic  processes. 

This  scientific  consortium  identified  the  Philippine  Sea  as  a  reasonable  venue  for  this 
effort.  This  region  is  bounded  to  the  south  by  the  North  Equatorial  Current  (NEC)  and  to  the 
west  by  the  Kuroshio.  Mesoscale  structures  propagate  westward  into  the  basin  and  collide 
with  eddies  spun  off  from  the  Kuroshio,  creating  energetic  and  complicated  oceanography. 

The  primary  scientific  experiment  is  planned  for  2010-2011,  and  involves  a  consider¬ 
able  number  of  new  systems.  To  mitigate  the  risk  of  deploying  new  systems,  the  group 
conducted  a  pilot  study /engineering  test  in  the  same  operating  environment  in  March-May 
2009. 

The  geographical  region  of  operations  and  the  assets  deployed  are  shown  in  Eig.  57. 
These  include: 


1.  An  autonomous  transceiver  mooring  (Tl).  This  unit  was  one  of  six  under  refurbish¬ 
ment  for  the  2010-201 1  experiment,  and  was  deployed  at  the  design  position  for  the 
Tl  mooring. 

2.  A  distributed  vertical  line  array  (VEA).  In  2009  this  was  a  5000-m  partially  populated 
development  version  of  a  system  targeted  for  the  main  experiment.  This  was  located 
at  the  design  intersection  of  two  transects  between  several  autonomous  transceiver 
moorings  scheduled  for  the  2010-2011  experiment. 

3.  Ship-suspended  stationary  operations  along  the  path  between  Tl  and  the  VEA.  Sta¬ 
tionary  operations  listed  herein  occurred  at  ship  stops  SS45  and  SS107. 


The  MP200/TR1446  was  freighted  to  Kao-Hsiung,  Taiwan,  in  March  2009,  and  lifted 
aboard  the  R/V  Melville  for  use  during  the  ONR  2009  Philippine  Sea  engineering  test  / 
pilot  study. 
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Figure  57:  Geographical  context  for  the  Philippine  Sea  2009  exercise.  APL-UW  deployed 
the  MP200/TR1446  at  ship  stops  SS45  and  SS107. 
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Figure  58:  Dockside  mobilization,  Kao-Hsiung  harbor,  April  2009.  The  MP200/TR1446 
system  was  transported  to  Kao-hsiung  in  a  softtop  freight  container. 


The  APL-UW  multiport  system  was  deployed  four  times  during  the  2009  exercise  — 
three  times  at  station  SS45  (no  transmissions  were  conducted  during  one  of  those  deploy¬ 
ments)  and  once  at  station  SS 107.  A  summary  of  transmissions  is  provided  in  Appendix  C. 

Although  the  2009  Philippine  Sea  exercise  was  an  engineering  test  of  equipment,  APL- 
UW  also  conducted  an  opportunistic  program  of  high  duty  cycle  full  power  transmissions 
to  the  Scripps  VLA  in  between  other  activities.  In  total,  the  program  netted  6.2  hours  of 
cumulative  transmission  time  out  of  9.5  hours  at  station  SS45,  and  32.2  hours  out  of  60 
hours  at  station  SS107.  A  comprehensive  analysis  of  those  transmissions  is  the  subject  of 
future  research. 


7.2  System  Limitations 

How  hard  can  the  system  be  driven?  Engineering  measurements  identified  several  con¬ 
straints. 
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Figure  59:  The  MP200/TR1446  and  the  HX554  on  the  fantail  of  the  R/V  Melville  prior 
to  departure  from  Kao-Hsiung  harbor.  The  SeaBattery  and  telemetry  bottle  are  not  yet 
attached  to  the  MP200/TR1446. 
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Figure  60:  Deploying  the  MP200/TR1446. 


Figure  61:  The  MP200/TR1446  descending  into  the  sea. 
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Figure  62:  Recovering  the  MP200/TR1446.  Two  hand  lines  and  two  lines  attached  to  air 
tuggers  were  used  to  control  the  system. 


Figure  63:  Attaching  the  monitor  hydrophone  and  its  bungy  mounting  frame  to  the  suspen¬ 
sion  cable  over  the  stern  of  the  R/V  Melville  during  deployment. 
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1 .  Krohn-Hite  Smoothing  Filter.  Both  sides  of  the  Krohn-Hite  filter  were  found  to  clip 
the  positive  side  of  an  input  waveform  above  about  8V;  the  negative  side  was  good  to 
—  lOV.  This  places  an  upper  limit  on  the  signal  file  amplitude  at  about  1500  quanta. 

2.  Instruments,  Inc.,  Amplifier.  The  Instruments,  Inc.,  L50  power  amplifier  was  mod¬ 
ified  at  APL-UW  to  convert  5  Vrms  to  1200  Vrms  on  the  voltage  source  setting 
1200/5.4.  This  is  the  setting  used  to  determine  the  conversion  scale  factor  of  258 
used  in  section  5.  Using  a  maximum  peak  value  of  1200  x  v/2  or  1697  V  and  a 
gain  of  258,  the  maximum  input  peak  voltage  is  6.58  V,  which  corresponds  to  a  peak 
signal  file  amplitude  of  1348  quanta.  This  places  a  second,  and  more  restrictive,  limit 
on  the  input  signal  amplitude. 

3.  MP200.  George  Cavanagh  of  Massa  advised,  “The  transducer  should  be  used  with 
the  equalizer  provided  to  ensure  you  do  not  overdrive  the  unit  at  the  high  Q  peaks 
in  the  response.  The  Transmit  Voltage  Response  (TVR)  between  the  peaks  appears 
...  to  be  about  130  dB.  A  Sound  Pressure  Level  (SPL)  of  195  dB  should  be  obtained 
at  about  1800  V  rms.”  [pers.  comm.,  10  October  2008]  This  corresponds  to  a  peak 
voltage  of  2550  V.  Given  the  step-up  autotransformer/tuner  and  the  power  amplifier 
gain,  this  corresponds  to  the  most  restrictive  limit  on  input  signal  amplitude  of  101 1 
quanta. 

7.3  Lessons  Learned 

Several  surprises  were  encountered  with  the  multiport  system  during  the  Philippine  Sea 
exercise. 


•  The  difference  between  the  theoretical  and  measured  scaling  between  the  monitor 
hydrophone  and  the  source  level  was  larger  than  expected.  The  theoretical  value 
was  given  (section  5.4.3)  as  30.4  dB,  but  the  value  measured  in  Lake  Washington  in 
January  2009  was  35.3. 

•  The  radiated  source  levels  were  generally  around  188.0  dB  for  the  first  few  transmis¬ 
sions  at  station  SS45  (see  Appendix  C),  much  less  than  predicted.  The  reason  for  the 
discrepancy  remains  unknown. 

•  The  transmission  levels  dropped  by  about  half  a  dB  after  the  first  few  transmissions 
at  station  SS45,  and  remained  there  the  rest  of  the  station.  The  reason  for  this  change 
also  remains  unknown. 

•  The  seawater  intrusion  into  the  tuner  housing  caused  irreparable  damage  to  the  tuner 
at  station  SS45.  The  work-around  was  to  wire  past  the  tuner,  and  drive  the  multiport 
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Figure  64:  Pulse-compressed  result,  transmission  file  1239933870.  sam,  as  measured 
on  the  monitor  hydrophone. 


transducer  directly  from  the  end  of  the  EOM  cable.  This  does  not  appear  to  have 
harmed  the  hardware.  The  loss  of  the  1:2  step-up  tuner  should  have  contributed  to  a 
6  dB  loss  in  radiated  source  level.  Measurements  based  on  the  monitor  hydrophone 
indicate  that  the  levels  were  consistently  about  182.5,  or  about  5  dB  lower.  Again, 
the  discrepancy  (here,  of  1  dB)  is  not  understood. 

•  The  pulse-compressed  waveforms  from  either  the  monitor  hydrophone  (Fig.  64)  or 
from  the  vertical  line  array  (not  shown)  suggest  that  the  equalizer  was  not  optimally 
tuned. 

As  an  example,  consider  the  pulse-compressed  output  (Fig.  64)  from  the  first  full 
power  transmission  from  SS45.  The  data  are  from  file  1239933870.  sam,  start 
time  UTC  107 :02:04: 1 1 .  The  data  were  captured  on  the  monitor  hydrophone,  located 
vertically  21  m  from  the  acoustic  center  of  the  MP200/TR1446.  The  propagation 
time  over  21  m  is  approximately  14  ms:  the  pulse  appears  to  arrive  at  about  16  ms, 
so  this  is  the  direct  path.  There  was  no  detectable  pulse  in  the  record  at  the  time 
expected  for  a  surface-reflected  path.  (The  depth  was  approximately  1000  m.)  The 
pulse  here  has  considerable  trailing  sidelobes,  similar  to  but  not  as  prominent  as  those 
in  the  uncompressed  pulse  example  (Fig.  54). 

Suboptimal  tuning  can  also  be  inferred  from  Fig.  56,  which  shows  results  from  the 
2009  Lake  Washington  test.  In  this  figure,  the  nulls  in  the  measured  spectrum  do 
not  align  with  those  in  the  theoretical  spectrum.  The  “ringing”  was  not  detected  in 
the  2009  Lake  Washington  data  because  any  “tail”  behind  the  direct  path  pulse  was 
obscured  by  the  arrival  of  the  surface-reflected  pulse  (see  Lig.  55.) 

While  there  should  exist  some  variation  in  the  radiated  signal  spectrum  from  broad¬ 
side  to  endfire,  this  variation  remains  unlikely  to  be  the  primary  reason  for  the  trailing 
sidelobes,  and  hence  further  refinements  in  the  equalizer  should  be  investigated. 
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It  is  no  surprise  that  the  performance  predictions  made  with  the  system  model  were 
not  completely  accurate.  Many  simplifying  approximations  were  made  to  establish  the 
system  model.  Nevertheless,  the  discrepancies  were  not  severe  considering  no  full  or  even 
moderate  power  calibration  was  conducted  with  the  fully  assembled  system.  Without  such  a 
calibration,  the  effects  on  acoustic  parameters  such  as  resonance  frequencies  and  broadside 
radiated  source  level  of  various  system-level  mechanical  issues,  such  as  the  tuner  mass  and 
SeaBattery  mass  mounted  at  the  transducer  equator,  remain  unknown. 

In  spite  of  all  the  engineering  approximations,  the  multiport  system  provided  many 
hours  of  usable  acoustic  data  during  the  2009  Philippine  Sea  exercise.  Further  refinements 
are  anticipated  for  future  experiments. 
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A  2006  Lake  Washington  Test  Results 


Twenty-one  test  measurements  were  conducted  during  the  Lake  Washington  test  of  2006. 
A  summary  of  each  collection  is  provided  in  Table  14. 

Results  of  the  admittance  calculations  for  each  transmission  are  shown  only  for  the 
frequencies  in  [100, 400]  Hz.  The  admittance  estimates  outside  this  range  were  too  noisy  to 
be  useful.  Occasional  frequency  points  are  annotated  on  the  figures  themselves.  The  name 
of  the  data  file  is  also  provided  in  the  figures. 


file 

start  time  (UTC) 

depth 

mptestOl.sam 

Apr  11  23:26:52  2006 

Depth  150  ft 

mptest02.sam 

Apr  11  23:34:49  2006 

Depth  150  ft 

mptestOS.sam 

Apr  11  23:41:59  2006 

Depth  150  ft 

mptest04.sam 

Apr  11  23:46:54  2006 

Depth  150  ft 

mptestOS.sam 

Apr  11  23:58:51  2006 

Depth  150  ft 

mptest06.sam 

Apr  12  00:05:16  2006 

Depth  150  ft 

mptestOV.sam 

Apr  12  18:27:00  2006 

Shallow 

mptestOS.sam 

Apr  12  18:34:31  2006 

Shallow 

mptest09.sam 

Apr  12  18:52:14  2006 

Shallow,  Law  was  103,  not  used  here 

mptestlO.sam 

Apr  12  18:59:51  2006 

Shallow,  Law  was  103,  not  used  here 

mptestl  l.sam 

Apr  12  19:06:05  2006 

Shallow 

mptestl2.sam 

Apr  12  20:14:12  2006 

Depth  100  ft 

mptestl  3. sam 

Apr  12  20:29:34  2006 

Depth  100  ft 

mptestl  4. sam 

Apr  12  20:35:14  2006 

Depth  100  ft 

mptestl  5. sam 

Apr  12  20:44:12  2006 

Depth  100  ft 

mptestl  6. sam 

Apr  12  20:48:33  2006 

Depth  100  ft 

mptestl  7. sam 

Apr  12  22:37:57  2006 

Dockside,  in  air:  garbage 

mptestl  8. sam 

Apr  12  22:43:55  2006 

Dockside,  in  air:  garbage 

mptestl  9. sam 

Apr  12  22:49:59  2006 

Dockside,  in  air:  garbage 

mptest20.sam 

Apr  12  23:09:36  2006 

Dockside,  1  m  deep  :  not  used  here 

mptest2  l.sam 

Apr  12  23:19:40  2006 

Dockside,  1  m  deep  :  not  used  here 

Table  14:  Summary  of  transmission  collections,  Lake  Washington  2006  test.  Most  trans¬ 
missions  used  an  RMS  drive  voltage  of  about  100  V. 
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B  System  SPICE  Netlist 

The  following  netlist  is  the  version  generated  automatically  by  the  simulator. 


automatic  netlist,  following  TechMemo:  2009.08.01 

*  D/A  subcircuit 

. subckt  DAconverter  1234 
rin  1  2  le+06 
eout  3412  0.00488 
.  ends 

■j'r'j'rVrVr'j'r'j'ryrVr'j'r'j'rVr^'j'r'j'ryr'j'r'j'r^'j'rVryr'j'r^^'j'ryr'j'r'j'r^'j'r 

*  Krohn-Hite  filter  stand-in: 

. subckt  KrohnHiteFilter  1234 

rKH  131 

cKH  3  2  3.18e-05 

.  ends 

*  ★VrVr'j'r'j'rVr'7'r'jir'jirVrVr'j'r'jicVr'jir'jirVrVcyrVrVrVrVr'j'rVr'jir'jirVr'j'r 


*  L50  power  amp: 

. subckt  PowerAmpLSO  1234 
rin  1  2  le+06 
el  3  4  1  2  258 
.  ends 


*  cable: 

. subckt  EOMcable  1234 
rl  1  5  3.15 
r2  5  3  3.15 
cl  5  2  1.05e-05 
.  ends 


*  tuner/autotransf ormer : 

. subckt  CoiltronTuner  123 

Ipl  3  5  0.15 

rpl  511 

lp2  1  6  0.15 

rp2  621 

kcoupling  Ipl  lp2  0.98 
.  ends 


•?'r'?k'VryrVr'?'rVryr'?'r'?'r':^yr'?'r'?iryr'?'r'?'r':^VrVryr'?'r'?irVr'?'rVr'?'r'?k'Vr'?'r 

*  transducer  MP  model,  S  =  8.12e4 
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. subckt  MPtransducer  12345 

cO  1  2  1.25e-06 

fl  1  2  vsensea  -1.23e-05 

el  6  5  1  2  1.23e-05 

vsensea  6  7  dc  0 

do  7  8  712 

llo  8  9  7.92e-10 

rlo  9  3  1.93e-08 

chi  7  10  181 

Ihi  10  11  1.37e-09 

rhi  11  4  1.06e-07 

.  ends 


*  circuit: 

vin  1  0  ac  1.00  dc  0.0 
xl  1  0  2  0  DAconverter 
x2  2  0  3  0  KrohnHiteFilter 
x3  3  0  4  0  PowerAmpLSO 
x4  4  0  5  0  EOMcable 
x5  5  0  6  CoiltronTuner 
x6  6  0  7  8  0  MPtransducer 
vsenselo  7  0  dc  0 
vsensehi  8  0  dc  0 


.ac  lin  12499  0.2  2499.8 
.options  NOPAGE 
.print  ac  i (vsensehi) 

.  end 
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C  2009  Philippine  Sea  Transmissions 

C.l  Background 

Signals  were  transmitted  during  the  2009  ONR  Philippine  Sea  Pilot  Study /Engineering  Test 
at  two  “ship  stops,”  SS45  and  SS107.  The  total  number  of  multiport  system  transmissions 
at  stations  SS45  and  SS107  were  19  and  214,  respeetively.  The  respeetive  total  transmission 
times  were  6.2  hours  spanning  9.5  hours  and  32.3  hours  spanning  60  hours. 

A  transmission  file  was  generated  for  eaeh  transmission.  The  start  time  is  given  in  the 
form  YYY:HH:MM:SS  where  YYY  is  the  yearday  with  January  1  being  001,  all  in  UTC. 
The  times  listed  below  are  the  aetual  times  of  the  first  seeond  of  transmission,  as  deeoded 
from  the  IRIG  ehannel.  The  marktime  is  20  seeonds  later,  and  may  be  different  by  1  seeond 
from  the  target  marktime  reeorded  in  the  file  header.  (See  below  for  explanation.)  The  ramp 
in  all  eases  is  20  seeonds. 


C.2  Station  SS45 

There  were  three  deployments  of  the  MP200/TR1446  at  SS45: 

1.  From  4/17  00:00Z  to  4/18  00:00Z  (day  107).  Seheduled  transmissions  were  inter¬ 
rupted  several  times  for  heat  build-up  issues  in  the  seience  van.  Transmissions  were 
aborted  due  to  shorting  in  the  underwater  paekage.  This  was  later  determined  to 
be  due  to  seawater  in  the  transformer  housing  (oil  expansion  system  tubing  diseon- 
neeted,  possibly  in  going  through  the  air-sea  interfaee  in  strong  eurrents.  ) 

2.  From  4/20  02:00Z  to  4/20  08:00Z  (day  110).  The  seeond  transmission  was  aborted 
due  to  shorting  in  the  underwater  paekage.  This  was  later  determined  to  be  a  leak  in 
the  splice:  there  was  also  possible  damage  to  the  autotransformer. 

3.  4/24  (day  114).  The  MP200/TR1446  was  deployed,  and  stopped  at  about  200  m  for 
a  check,  but  it  was  determined  that  the  underwater  package  was  shorted.  It  was  re¬ 
covered,  and  the  autotransformer  was  found  to  be  shorted.  No  transmissions  resulted 
from  this  deployment. 

The  signal  file  in  all  cases  was  m2 8 4  .EQ.wav,  with  parameters  given  in  Table  15. 
Transmissions  are  listed  in  Table  16. 

The  system  time  was  observed  to  be  drifting  relative  to  the  local  GPS  and  ship  NTP 
time  servers.  The  reason  for  the  drift  is  unknown  as  it  was  not  observed  in  the  laboratory 
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carrier 

284.166666  Hz 

cycles/digit 

2 

digit  length 

7.04  ms 

sequence  length 

1023 

law 

2033 

filtering 

eqOl.nc 

Table  15:  Signal  parameters,  2009  Philippine  Sea  experiment. 


at  APL-UW.  The  two  time  servers  did  not  differ  by  more  than  100  ms,  generally  much 
less.  The  effect  of  this  drift  is  that  the  actual  transmission  time  may  be  the  second  adjacent 
to  the  desired  start  time:  the  transmission  program  (which  was  reading  the  system  clock) 
triggered  on  the  desired  second  as  read  from  the  system  clock,  but  the  hardware  (which  was 
synchronized  by  the  local  GPS)  was  on  a  different  second,  i.e.,  based  on  cruise  hindsight, 
the  GPS  second  was  in  fact  the  desired  second,  and  the  time  read  from  the  system  clock  was 
incorrect.  At  no  time  during  occasional  checks  and  subsequent  testing  was  it  ever  observed 
that  the  GPS  lost  lock. 

Table  16:  Multiport  transmissions  for  station  SS45  in  the 

2009  Philippine  Sea  experiment.  Duration  includes  ramp. 


start  time 

duration 

source 

file  name 

[UTC] 

[s] 

level 

comments 

1239933870.sam 

107:02:04:11 

3046.2 

188.1 

off  by  1  second 

1239937038.sam 

107:02:56:58 

253.3 

188.0 

1239937470.sam 

107:03:04:10 

647.5 

188.1 

Aborted 

1239941070.sam 

107:04:04:10 

3047.2 

188.2 

1239944238.sam 

107:04:56:58 

253.3 

187.8 

123995 1438.sam 

107:06:56:58 

252.0 

190.9 

Sample  rate  error 

1239951870.sam 

107:07:04:10 

3033.3 

190.9 

Sample 

rate  error 

1239955038.sam 

107:07:56:59 

252.5 

187.7 

off  by  1 

second 

1239955470.sam 

107:08:04:11 

3046.5 

187.8 

off  by  1 

second 

1239958638.sam 

107:08:56:59 

252.5 

187.6 

off  by  1 

second 

1239959070.sam 

107:09:04:11 

3046.2 

187.7 

off  by  1 

second 

1239962238.sam 

107:09:56:59 

253.2 

187.5 

off  by  1 

second 

1239962670.sam 

107:10:04:10 

3047.2 

187.8 

1239965838.sam 

107:10:56:59 

253.2 

187.7 

off  by  1 

second 

1239966270.sam 

107:11:04:10 

1381.0 

187.3 

aborted 

1240203438.sam 

110:04:56:58 

253.3 

187.5 

1240203870.sam 

110:05:04:11 

3046.2 

187.6 

off  by  1 

second 

continued  on  next  page 

C2  TR  0902 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


file  name 

start  time 
[UTC] 

duration 

[s] 

source 

level 

comments 

1240207038.sam 

1240207470.sam 

110:05:56:59 

110:06:04:11 

252.5 

3020.8 

187.4 

187.5 

off  by  1  second 

off  by  1  second,  aborted 

C.3  Station  SS107 

The  MP200/TR1446  was  deployed  at  station  SS107  from  about  4/27  01:30Z  (day  117)  to 
about  4/29  16:00Z  (day  119).  The  autotransformer  was  electrically  disconnected  in  the 
underwater  package,  and  the  power  amplifier  signal  was  essentially  driven  directly  into  the 
MP200/TR1446. 

Initial  tests  at  depth  indicated  that  the  multiport  itself  would  only  require  about  1.6  A 
rms  at  the  maximum  drive  voltage,  so  we  felt  more  confident  that  heat  build-up  would  not 
be  a  problem.  Therefore,  we  scheduled  the  source  using  the  file  m2  8  4  .  EQ  .  wav  on  the 
normal  long/short  schedule.  This  worked  for  about  2  hours,  but  then  we  experienced  two 
failures.  The  first  was  a  change  in  cable  voltage  and  current,  along  with  a  corresponding 
increase  in  radiated  power.  This  was  likely  due  to  a  shift  in  D/A  trigger  clock  frequency. 
Both  the  time  board  and  GPS  were  replaced  by  backup  units.  The  other  problem  was 
sporadic  fault  conditions  on  the  power  amplifier.  This  problem  might  last  only  one  period, 
or  many  seconds.  The  only  known  corrective  procedure  at  the  time  was  to  give  the  cabinet 
a  shake  and  hope  for  the  best.  (This  worked  3  times  out  of  3.) 

Occasionally  a  transmission  was  aborted  manually  when  a  fault  came  on.  These  are 
noted  below.  Nevertheless,  based  on  the  notion  that  these  faults  were  exacerbated  by  heat 
build-up,  the  transmission  schedule  was  changed  to  one  closer  to  a  50%  duty  cycle.  To 
avoid  complete  transmission  down  time  while  the  new  schedule  was  being  recalculated 
and  verified,  several  10-minute  transmissions  were  launched  manually.  These  are  unlikely 
to  be  synchronized  to  the  original  m-sequence  schedule,  and  are  labeled  “disparate”  in 
Table  17. 

This  new  automatic  schedule,  utilizing  35  transmission  minutes  per  hour,  began  around 
10:00Z  on  day  1 17  (4/27).  This  schedule  ran  10  min  on  5  off,  10  on  5  off,  15  on  15  off  per 
hour.  This  schedule  was  used  throughout  the  night. 

At  approximately  01:00Z  on  118  (4/28)  the  schedule  was  modified  to  add  150  seconds 
of  transmission  data  at  the  top  of  the  hour.  This  was  primarily  to  get  a  burst  of  measure¬ 
ments  in  the  otherwise  long  gap  starting  about  48  minutes  after  the  hour. 

Amplifier  glitches  were  occasionally  observed  by  watchstanders:  these  lasted  from  one 
period  to  several  seconds  or  more.  The  files  have  not  been  cleansed  of  any  consequent 
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waveform  interruptions,  and  the  locations  of  such  anomalies  remain  to  be  determined  from 
the  files  themselves. 


Table  17:  Multiport  transmissions  for  station  SS107  in  the 
2009  Philippine  Sea  experiment.  Duration  includes  ramp. 


file  name 

start  time 
[UTC] 

duration 

[s] 

source 

level 

comments 

1240804596.sam 

117:03:56:16 

253.3 

182.8 

1240805028.sam 

117:04:03:28 

3047.3 

182.9 

1240808196.sam 

117:04:56:16 

253.3 

182.7 

1240808628. sam 

117:05:03:28 

3047.3 

182.9 

124081 1796.sam 

117:05:56:16 

251.2 

185.8 

samplerate  error 

12408 12228.  sam 

117:06:03:28 

3023.6 

185.8 

samplerate  error 

1240815396.sam 

117:06:56:16 

251.2 

185.6 

samplerate  error 

12408 15828.sam 

117:07:03:28 

765.8 

185.6 

samplerate  error, 
aborted? 

ss  107-disparateA.sam 

117:08:22:43 

621.3 

N/A 

manual  launch 

ssl07-disparateB.sam 

117:08:34:02 

621.3 

N/A 

manual  launch 

1240822596. sam 

117:08:56:16 

253.3 

182.7 

1240823028. sam 

117:09:03:28 

1271.8 

182.5 

aborted 

ssl07-disparateC.sam 

117:09:27:27 

621.3 

N/A 

manual  launch 

1240826628. sam 

117:10:03:28 

621.3 

182.7 

1240827528.sam 

117:10:18:28 

621.3 

182.7 

1240828428. sam 

117:10:33:28 

921.3 

182.7 

1240830228. sam 

117:11:03:28 

621.3 

182.7 

1240831 128.sam 

117:11:18:28 

621.3 

182.7 

1240832028. sam 

117:11:33:28 

921.3 

182.7 

1240833828. sam 

117:12:03:28 

621.3 

182.7 

1240834728.sam 

117:12:18:28 

621.3 

182.7 

1240835628. sam 

117:12:33:28 

921.3 

182.7 

1240837428.sam 

117:13:03:28 

621.3 

182.7 

1240838328. sam 

117:13:18:28 

621.3 

182.7 

1240839228.sam 

117:13:33:28 

921.3 

182.8 

1240841028.sam 

117:14:03:28 

621.3 

182.6 

1240841928.sam 

117:14:18:28 

621.3 

182.7 

1240842828.sam 

117:14:33:28 

921.3 

182.7 

1240844628.sam 

117:15:03:28 

621.3 

182.7 

1240845528.sam 

117:15:18:28 

621.3 

182.7 

1240846428.sam 

117:15:33:28 

921.3 

182.7 

1240848228.sam 

117:16:03:28 

621.3 

182.7 
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file  name 

start  time 
[UTC] 

duration 

[s] 

source 

level 

comments 

1240959828.sam 

118:23:03:28 

1240960728.sam 

118:23:18:28 

mm 

1240961 628.  sam 

118:23:33:28 

921.0 

182.4 

1 240962996.  sam 

118:23:56:16 

171.0 

182.3 

1240963428.sam 

119:00:03:28 

621.0 

182.4 

1240964328.sam 

119:00:18:28 

621.0 

182.5 

1240965228.  sam 

119:00:33:28 

921.0 

182.5 

1240966596. sam 

119:00:56:16 

171.0 

182.4 

1240967028.sam 

119:01:03:28 

621.0 

182.5 

1240967928.sam 

119:01:18:28 

621.0 

181.7 

1240968828. sam 

119:01:33:28 

921.0 

182.4 

1240970628.sam 

119:02:03:28 

621.0 

182.4 

1240971528.sam 

119:02:18:28 

241.0 

181.7 

Aborted 

1240972428.sam 

119:02:33:28 

921.0 

182.4 

1240973796.sam 

119:02:56:16 

171.0 

182.3 

1240974228.sam 

119:03:03:28 

621.0 

182.5 

1240975128.sam 

119:03:18:28 

621.0 

182.4 

1240976028.sam 

119:03:33:28 

921.0 

182.5 

1240977396.sam 

119:03:56:16 

171.0 

182.2 

1240977828.sam 

119:04:03:28 

621.0 

182.8 

1240978728.sam 

119:04:18:28 

621.0 

182.4 

1240979628.sam 

119:04:33:28 

921.0 

182.4 

1240980996.sam 

119:04:56:16 

171.0 

182.2 

1240981428.sam 

119:05:03:28 

621.0 

182.2 

1240982328. sam 

119:05:18:28 

621.0 

182.2 

1240983228.sam 

119:05:33:28 

921.0 

182.4 

1240984596.sam 

119:05:56:16 

171.0 

182.2 

1240985028.sam 

119:06:03:28 

621.0 

185.7 

1240985928.sam 

119:06:18:28 

621.0 

183.0 

1240986828. sam 

119:06:33:28 

921.0 

187.9 

1240988196.sam 

119:06:56:16 

171.0 

182.4 

1240988628. sam 

119:07:03:28 

621.0 

182.2 

1240989528.sam 

119:07:18:28 

621.0 

182.2 

1240990428.sam 

119:07:33:28 

921.0 

182.4 

1 24099 1796.sam 

119:07:56:16 

171.0 

182.1 

1240992228. sam 

119:08:03:28 

621.0 

182.3 

1240993128.sam 

119:08:18:28 

621.0 

182.4 

1240994028.sam 

119:08:33:28 

921.0 

182.3 
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file  name 

start  time 
[UTC] 

duration 

[s] 

source 

level 

comments 

1240995396.sam 

119:08:56:16 

171.0 

1240995828.sam 

119:09:03:28 

621.0 

mm 

1240996728.sam 

119:09:18:28 

621.0 

1240997628.sam 

119:09:33:28 

921.0 

182.3 

1240998996.sam 

119:09:56:16 

171.0 

182.1 

1240999428.sam 

119:10:03:28 

621.0 

182.3 

1241000328.sam 

119:10:18:28 

621.0 

182.3 

1241001228.sam 

119:10:33:28 

921.0 

182.4 

1241002596.sam 

119:10:56:16 

171.0 

182.2 

1241003028.sam 

119:11:03:28 

621.0 

182.3 

1241003928.sam 

119:11:18:28 

621.0 

182.2 

1241004828.sam 

119:11:33:28 

921.0 

182.2 

1241006196.sam 

119:11:56:16 

171.0 

182.0 

1241006628.sam 

119:12:03:28 

621.0 

182.2 

1241007528.sam 

119:12:18:28 

621.0 

182.2 

1241008428.sam 

119:12:33:28 

921.0 

182.3 

1241009796.sam 

119:12:56:16 

171.0 

182.3 

1241010228.sam 

119:13:03:28 

621.0 

182.2 

124101 1128.sam 

119:13:18:28 

621.0 

182.2 

1241012028.sam 

119:13:33:28 

921.0 

182.2 

1241013396.sam 

119:13:56:16 

171.0 

182.1 

1241013828.sam 

119:14:03:28 

621.0 

182.4 

1241014728.sam 

119:14:18:28 

621.0 

182.3 

1241015628.sam 

119:14:33:28 

921.0 

184.3 

1241016996.sam 

119:14:56:16 

171.0 

182.1 

1241017428.sam 

119:15:03:28 

621.0 

182.1 

1241018328.sam 

119:15:18:28 

621.0 

182.3 

1241019228.sam 

119:15:33:28 

921.0 

183.2 

1241020596.sam 

119:15:56:16 

171.0 

182.4 

1241021028.sam 

119:16:03:28 

87.8 

165.4 
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